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ABSTRACT 


Two  topics  dealing  with  the  reflection  of  ultrasound  bursts  from 
random  media  are  discussed  In  this  work. 

In  chapter  2  we  develop  a  general  formulation  of  the  echo  received 
from  a  random  scatterer  ensemble  Illuminated  by  a  short  electromagnetic 
or  sonic  signal.  We  show  theoretically  that  a  gradient  In  either  scat¬ 
terer  concentration  or  In  the  field  function  of  the  transmitter/receiver 
will  return  an  echo  which  Is  partially  spatially  coherent  (l.e.  specular). 
Furthermore  we  show  that  from  the  degree  of  coherency,  l.e.  from  the  ratio 
of  the  random  part  to  the  nonrandom  part  of  the  reflected  signal,  the 
scatterer  concentration  and  scattering  cross-section  can  be  calculated. 
We  also  show  experimentally  that  scattering  concentration  gradient  creates 
a  coherent  reflection  from  whose  degree  of  coherency  the  scattering  con¬ 
centration  can  be  estimated. 

In  chapter  3  and  A  we  deal  with  signal  processing  techniques  for  the 
reduction  of  clutter  noise  created  by  random  media  composed  of  a  high 
concentration  of  point  scatterers.  The  problem  Is  to  enhance  a  target 
embedded  In  a  random  medium  when  the  clutter  noise  variance  and  the 
target  echo  location  and  amplitude  are  unknown.  In  chapter  3  we  analyze 
theoretically  a  technique  that  was  first  proposed  by  Newhouse  et  al  [3] 
the  so  called  Minimization  algorithm.  In  this  technique  we  split  the 
received  signal  spectrum  Into  n  frequency  windows.  The  minimum  of  the 
squared  signals  at  each  range  delay  Is  then  chosen.  The  calculated  SNRE 
(slgnal-to-nolse-ratlo  enhancement)  of  this  technique  agrees  well  with 
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experiments  performed  by  N.M.  Bilgulay  [4] .  We  calculate  also  the 
Receiver  Operating  Characteristics  (ROC)  and  find  out  that  the  Improved 
SNRE  of  Minimization  Is  negated  by  loss  in  detection  properties. 

In  chapter  4  we  describe  a  method  for  clutter  reduction  that  also 
uses  the  split  spectrum  principle.  In  this  technique  we  construct  the 
optimum  receiver  for  each  range  delay  from  the  n  outputs  of  the  frequency 
windows.  We  calculate  the  ROC  for  this  technique  and  find  it  to  be 
superior  to  minimization  for  most  cases.  We  show  experimentally  the 
effectiveness  of  this  algorithm  In  clutter  reduction. 

In  chapter  5  we  summarize  the  results  of  this  research  work  and  show 
the  flexibility  and  versatility  of  the  techniques  by  introducing  Improve¬ 
ments  and  applications  to  the  described  methods. 
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CHAPTER  1 


INTRODUCTION 

CLUTTER  PROBLEMS  IN  IMAGING  SYSTEMS 

One  of  the  most  Important  limitations  of  Ultrasonic  Imaging  Systems 
using  pulse  echo  techniques  Is  Imposed  by  clutter  noise  at  the  system 
output.  This  clutter  noise  Is  usually  caused  by  relatively  small,  highly 
dense,  randomly  positioned  scatterers.  In  one  type  of  clutter  problem 
the  target  to  be  detected  can  be  modelled  as  a  strong  reflector.  Even 
though  the  target  echo  Is  significantly  larger  than  each  of  the  Individual 
small  random  scatterers  around  It,  the  target  Is  sometimes  difficult  to 
detect  due  to  the  high  density  of  the  Interfering  scatterers.  Such  prob¬ 
lems  exist  in  almost  every  field  associated  with  Imaging.  In  Ultrasound 
Nondestructive  Testing  the  large  grain  boundary  echoes  makes  it  dlf  ~  \lt 
to  detect  relatively  large  flaws  [1,2, 3, 4].  In  medical  Imaging  the  fine 
tissue  microstructure  echoes  sometimes  makes  It  hard  to  outline  the  organ 
boundary  [5,6],  and  In  Radar  Systems,  rain  drops,  chaff  or  sea  clutter 
[7,8,9,10,11,12,13,14]  are  often  strong  enough  to  mask  the  target.  Time 
averaging  or  correlation  techniques,  which  reduce  random  thermal  noise 
significantly,  are  not  suitable  for  reducing  coherent  noise  resulting 
from  echoes  due  to  the  stationary  Interfering  scatterers. 

In  another  type  of  clutter  problem  the  targets  are  composed  of  a 
high  density  of  point  scatterers.  The  reflected  signal  Is  accompanied 
by  a  speckle  pattern  (spatial  brightness  fluctuations)  which  is  due  to 
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the  random  Interference  pattern  of  the  reflectors  composing  the  target. 
This  type  of  problem  appear  often  In  coherent  optical  Imaging  [15]  when 
the  object  (that  can  be  modelled  by  a  high  density  of  small  random  re¬ 
flectors)  Is  Imaged  by  a  highly  coherent  Illumination.  In  Ultrasound 
Medical  Imaging  the  echo  from  the  tissue  microstructure  which  can  be 
modelled  In  a  similar  way.  Is  an  Important  part  of  the  Image  and  helps  to 
dlferentlate  and  characterize  the  different  organs  [16,17,18,5,6,19-21, 
22-25]  . 

In  the  text  we  refer  to  the  first  problem  in  which  we  wish  to  sup¬ 
press  the  clutter  echo  and  enhance  a  large  target  echo  as  a  problem  In 
“clutter  reduction”,  and  we  refer  to  the  second  problem  in  which  we  wish 
to  reduce  the  speckle  pattern  of  targets  composed  of  a  large  sum  of  point 
reflectors  as  a  problem  In  "speckle  reduction”. 

Some  conventional  techniques  to  reduce  clutter  and  speckle  are 
described  below. 

GENERAL  LITERATURE  REVIEW 
Scattering  from  Random  Media 

Ultrasonic  signal  processing  techniques  have  developed  dramatically 
In  the  last  decade  mainly  due  to  the  Introduction  of  low  cost  powerful 
computers.  But  before  reviewing  the  recent  literature  on  modern  signal 
processing  techniques  In  Ultrasound  we  should  pay  some  attention  to  the 
basic  understanding  of  the  scattering  from  random  media  as  the  modelling 
of  random  media  Is  the  basis  on  which  most  of  the  signal  processing 
techniques  base  their  approach. 

The  theory  covering  the  propagation  of  waves  In  random  media  Is 
vast.  In  [26,27,28]  one  can  find  the  description  of  the  basic  acoustical 
models  for  wave  propagation  and  In  [29,30]  a  discussion  on  some  of  the 
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theories  when  applied  to  biological  tissues  Is  given. 

The  scattering  from  bodies  with  simple  geometry  can  be  calculated  by 
methods  Introduced  In  [26,27,31].  In  the  limit,  for  reflectors  which  are 
much  smaller  than  the  wave  length  the  solutions  are  explicit.  These 
solutions  can  be  extended  to  complex  elements  by  assuming  a  combination 
of  several  Individual  centers. 

The  scattering  Is  described  through  the  definition  of  scattering 
cross  section.  The  total  cross  section  Is  defined  as  [26,32] 

a  -  I  1.1 

I 

where  I  Is  the  Incident  wave  and  P  Is  the  scattered  signal.  Of  more 
practical  Importance  Is  the  definition  of  differential  scattering  cross 
section  which  la  related  to  o  by 

a  -  /  a’dA  1.2 

4ir 

where  A  Is  the  solid  angle.  In  general  a'>a'(A).  For  the  Important 
case  of  a  small  sphere  (smaller  than  the  wave  length  [26]) 

a'  -  i  [ll]^  a6[^ +  8  cos  9]2  If  X  »  a 

9  X  K  2p8+P 

where 

X  *  wavelength 
a  -  radius  of  the  scatterer 

9  *  angle  between  the  Incident  wave  and  the  direction  of  scattering. 

Pg  ■  density  of  the  sphere, 
p  ••  density  of  the  medium. 

Kg  >  compressablllty  of  the  sphere. 

K  compressablllty  of  the  medium. 


< 
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For  practical  cases  Isotropic  scattering  can  be  assumed  only  If 
either  Ps^p  or  If  the  range  of  angles  In  eq.  3  Is  small  enough  and  o  Is 
such  that  o-a'A^. 

The  structure  of  biological  tissue  Is  very  complex  and  It  exhibits 
nonhomogeneities  In  almost  any  scale.  So  the  assumption  of  discrete,  weak 
scattering  which  Is  often  used  Is  clearly  a  simplified  assumption  as  the 
scattering  bodies  range  from  much  smaller  than  the  wave  length  up  to  much 
larger  than  the  range  cell  [32] •  Theoretical  results  based  on  these 
assumptions  should  be  constantly  challenged  and  verified  to  be  considered 
reliable. 

The  statistical  properties  of  the  backscattered  clutter  echo  play  a 
maj^or  role  In  the  optimal  signal  processing  to  be  chosen  for  clutter  or 
speckle  reduction.  It  can  be  shown  that  If  the  scatterers  are  uniformly 
distributed  and  the  number  of  scatterers  In  the  range  cell  Is  high  enough, 
the  backscattered  echo  amplitude  assuming  the  plane  or  spherical  wave 
approximation,  can  be  considered  Gaussian  with  zero  mean  [33].  In  [34] 
the  average  power  backscattered  from  random  media  with  a  density  profile 
Is  given  for  continuous  wave  propagation.  In  this  publication  Slegert 
and  Goldstein  show  that  a  medium  which  exhibits  a  density  profile  gives 
rise  to  a  coherent  term  which  translates  to  an  echo  with  non-zero  mean. 
In  [35,36]  Glotov  calculates  the  average  backscattered  power  from  a  slab 
filled  with  nonuniform  size  scatterers  Illuminated  by  a  burst.  The  echo 
again  contains  a  coherent  term  which  resembles  specular  reflection.  If 
the  slab  Is  much  smaller  than  the  wave  length  the  reflection  Is  mostly 
coherent  and  the  random  component  becomes  Insignificant.  In  chapter  2 
we  will  develop  a  general  formulation  for  the  backscattered  echo  from 
random  media  on  which  our  signal  processing  approaches  will  be  based. 
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Signal  Processing  Techniques  for  Clutter  Reduction 

One  of  the  most  popular  techniques  used  In  radar  systems  to  overcome 
clutter  noise  Is  the  use  of  frequency  diversity  (agility).  In  this  tech¬ 
nique  the  Radar  system  possesses  a  center  frequency  which  Is  varied  between 
pulses.  This  Is  similar  to  Increasing  the  effective  bandwidth  of  the 
transmitted  signal  and  thus,  the  clutter  noise  power  which  Is  Inversely 
proportional  to  the  system  bandwidth  [7,8,14],  Is  reduced.  The  other 
Improvements  obtained  with  the  use  of  frequency  agile  Radars  namely, 
antijamming  capabilities,  range  Improvement  and  tracking  are  not  of  any 
Importance  in  biomedical  Ultrasound  and  non-destructive  testing.  However, 
techniques  used  for  processing  the  frequency  diverse  signals  will  be  the 
measure  against  which  the  new  techniques  (to  be  Introduced  In  chapters  3 
and  4)  will  be  compared. 

A  very  good  review  of  the  literature  on  frequency  agility  Is  given 
In  [4].  To  summarize,  we  only  mention  that  the  two  most  popular  proces¬ 
sing  schemes  are  composed  of  either  direct  summation  of  the  successive 
echoes  before  demodulation  (which  Is  equivalent  to  operation  with  a  wide 
bandwidth  or  transmitting  and  receiving  with  all  the  frequency  windows 
simultaneously)  or  summing  the  different  frequency  windows  after  demodula¬ 
tion  (envelope  detection).  This  Is  done  when  the  target  echo  phase  Is 
lost  between  pulses  (due  to  fast  fluctuation  of  the  target  position 
between  pulses). 

After  processing  the  frequency  diverse  signals  a  threshold  is  set  to 
enable  a  decision  whether  a  target  exists  In  a  prespecified  range  cell. 
In  the  technique  which  is  most  often  used  and  called  "CELL  AVERAGING"  we 
evaluate  the  clutter  noise  power  [9,37,38,10,11,12,13]  from  some  "test 
cells".  Then  assuming  that  the  same  clutter  noise  power  exists  In  the 
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range  cell  of  interest  a  CFAR  (constant  false  alarm  rate)  detector  can 
be  designed-  This  technique  provides  a  noise  riding  threshold  and  only 
cells  assumed  to  contain  targets  are  presented.  The  detector  design 
depends  on  the  statistics  of  the  clutter  echo.  Usually  the  assumption 
of  Gaussian  distribution  for  the  clutter  noise  amplitude  is  used  which 
is  true  only  if  the  number  of  scattering  centers  in  the  range  cell  is 
high  enough.  A  good  review  of  clutter  properties  (in  Radar  applications) 
Is  given  in  [39]. 

In  1981  "Split  Spectrum  Processing"  was  introduced  by  Newhouse  et  al 
[3,4]  for  the  Improvement  of  flaw  visibility  in  the  presence  of  grain 
noise.  In  this  technique  the  received  echo  is  split  into  several  fre¬ 
quency  windows  (possibly  overlapping)  and  the  minimum  of  the  squared 
windowed  signals  at  each  range  delay  is  then  chosen.  In  [3,4]  the  tech¬ 
nique  is  compared  experimentally  with  other  conventional  techniques  and 
proved  to  have  an  Improved  flaw  to  grain  noise  enhancement.  As  the 
'Optimal  Processing'  technique  introduced  in  Chapter  4  in  this  report  is 
partially  based  on  elements  in  this  algorithm  and  will  be  compared  with 
it,  we  analyze  the  minimization  algorithm  theoretically  in  chapter  3. 

Signal  Processing  Techniques  for  Speckle  Reduction  in  Ultrasonic  Imaging 

Because  the  Laser  speckle  phenomenon  is  so  similar  to  Ultrasound 
speckle  the  Laser  literature  was  the  main  source  for  techniques  for 
speckle  reduction  in  Ultrasound.  An  excellent  review  of  the  optical 
speckle  phenomenon  is  given  in  [13].  It  was  only  towards  the  middle  of 
the  last  decade  that  techniques  for  Ultrasound  speckle  reduction  were 
introduced  and  a  substantial  number  of  signal  processing  techniques 
appeared  in  the  open  literature  especially  in  the  last  few  years  [17,18,5, 
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20,22-24].  One  of  the  earliest  papers  [40]  uses  spatial  compound  scanning 
to  reduce  Ultrasound  speckle.  In  this  technique  the  region  of  Interest 
Is  Illuminated  from  different  directions  and  the  different  gray  levels 
are  then  compounded  to  create  an  Image.  This  technique  can  be  successful 
only  when  the  Interrogated  medium  can  be  considered  Isotropic  and  non 
refractive  and  the  region  of  Interest  Is  accessible  from  different  direc¬ 
tions.  Organs  such  as  the  breast  are  Ideal  for  such  compound  Imaging. 
In  1978  Burckhardt  argued  for  the  first  time  that  Ultrasound  speckle 
should  be  treated  In  a  similar  way  to  optical  speckle.  He  modelled  the 
tissue  as  a  collection  of  a  high  density  of  point  scatterers  and  he 
evaluated  the  signal  to  noise  ratio  as  the  mean  of  the  signal  envelope 
to  the  envelope  standard  deviation.  From  this  definition  the  signal  to 
noise  ratio  is  1.91.  He  later  Introduces  a  new  algorithm  "Compound  scan 
with  maximum  amplitude  writing".  In  this  technique  the  maximum  of  the 
different  echoes  associated  with  a  certain  point  In  space  Is  selected 
Instead  of  the  average.  He  claims  that  the  difference  between  averaging 
which  Is  the  optimum  processing  for  this  model  and  "Maximization"  Is 
negligible. 

In  1979  Abott  and  Thurstone  [19]  analyzed  the  Ultrasound  speckle 
phenomenon  based  again  on  the  concept  of  Laser  speckle.  The  paper  out¬ 
lines  the  ways  In  which  speckle  can  be  reduced  l.e  variations  In  time, 
space  or  frequency  parameters.  The  first  of  these,  time.  Involves  a 
phase  diffuser  which  Is  equivalent  to  generating  a  random  pattern  In 
phase  from  burst  to  burst  or  Image  to  Image.  This  technique  was  found 
to  be  useless  In  Ultrasound  as  It  resulted  In  Image  degradation  due  to 
deterioration  In  focusing  and  distortion  of  the  transmitted  wave  front. 
The  second  of  these,  spatial  variation.  Involves  either  illuminating  the 
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object  from  different  directions  or  simply  moving  the  transducer  in 
relation  to  the  object  sufficiently  to  change  the  phase  distribution 
across  the  transducer  aperture*  The  third  of  these  Is  equivalent  to 
frequency  agility  a  concept  that  was  discussed  earlier  In  connection 
with  target  detection  and  will  be  discussed  later  In  this  section  In 
connection  with  speckle  reduction. 

In  1983  Wagner  et  al  published  a  pair  of  papers  on  speckle.  In  the 
first  [16]  they  calculate  the  second  order  statistics  of  Ultrasound 
speckle  In  B  Images  as  a  function  of  transducer  dimensions  and  range. 
Their  results  seem  to  fit  reasonably  with  experiment.  In  the  second  paper 
[17]  they  use  the  results  obtained  In  the  first  paper  and  some  known 
results  from  statistical  decision  theory  to  obtain  decision  rules  for 
lesion  detection.  The  results  turn  out  to  be  spatial  averaging.  Several 
adjacent  cells  are  summed  to  create  a  new  cell  (or  a  decision  rule). 
Obviously  resolution  Is  sacrificed  In  this  procedure.  In  1983  [5]  Robinson 
and  Knight  checked  experimentally  the  performance  of  some  spatial  pulse 
echo  compound  scan  techniques.  The  performance  of  peak  detected,  minimum 
detected  and  averaged  reconstruction  of  point  targets  Is  compared.  They 
conclude  that  averaging  Is  capable  of  Increased  range  resolution  compared 
to  peak  detection  and  that  the  speckle  pattern  Is  smoother.  The  minimum 
detected  signal  has  an  Improved  resolution  of  point  targets  and  "shows 
promise  In  location  of  shadowing  structures”. 

Frequency  Compounding 

A  discussion  on  frequency  compounding  In  connection  with  Ultrasonic 
medical  Imaging  Is  given  In  [22].  It  Is  found  that  the  degree  of  speckle 
contrast  reduction  is  Inversely  proportional  to  the  bandwidth  of  the 
transmitted  acoustic  burst.  Also,  a  considerable  increase  In  signal  to 


noise  ratio  of  the  speckle  pattern  can  be  achieved.  In  another  paper 
[23]  the  received  spectrum  of  the  backscattered  signal  is  split  into  two 
overlapping  frequency  windows  and  the  respective  channels  are  then  summed 
after  envelope  detection.  This  technique  is  reported  to  Increase  the 
signal  to  noise  ratio  of  the  speckle  pattern  with  Improved  resolution. 
We  will  show  later  in  the  this  report  that  the  optimum  receiver  is  some¬ 
what  different  from  the  one  implemented  in  [23]  and  should  provide  better 
results. 

MOTIVATION  FOR  RESEARCH 

The  fact  that  current  techniques  are  not  effective  enough  in  the 
reduction  of  clutter  and  speckle  noise  in  Ultrasonic  imaging  systems  is 
the  driving  force  for  this  project.  Techniques  such  as  spatial  averaging 
(or  spatial  compounding)  are  difficult  to  Implement  due  to  the  relatively 
high  refraction  of  the  biological  tissue  and  man-made  materials,  the  high 
accuracy  required  from  the  scanning  system,  and  the  fact  that  many  of  the 
regions  of  Interest  are  not  accessible  from  different  directions.  Fre¬ 
quency  diversity  (or  frequency  averaging)  la  successful  when  it  is  easier 
to  transmit  narrow  frequency  windows  one  at  a  time  Instead  of  the  entire 
available  spectrum.  In  Ultrasound  however  it  is  possible  to  transmit 
extremely  wide  bandwidths. 

Extending  the  frequency  range  by  using  several  elements  with  differ¬ 
ent  center  frequency  raises  more  problems  than  it  solves.  First  the 
propagating  medium  la  highly  frequency  dependent  (l.e  frequency  attenua¬ 
tion  dependence)  and  secondly  the  scattering  from  the  microstructure 
rapidly  Increases  with  frequency.  The  beam  pattern  of  the  transducer  is 
also  sensitive  to  frequency.  Thus  simple  frequency  compounding  is  not 
useful. 


In  this  work  we  concentrate  our  effort  on  a  different  type  of  signal 
processing.  We  concentrate  on  clutter  and  speckle  reduction  through  post 
reception  A  scan  processing  only.  By  working  on  individual  A  scans  we 
eliminate  the  many  practical  problems  associated  with  compound  scanning 
and  make  the  problem  solely  a  signal  processing  one.  The  techniques  to 
be  Introduced  here  require  that  the  signals  should  be  processed  at  the  RF 
level.  Techniques  usually  described  in  the  literature  do  not  deal  with 
processing  at  the  RF  level,  but  as  we  see  later  the  hardship  of  working  at 
the  RF  level  has  its  rewards. 

SYNOPSIS  OF  THE  REMAINING  CHAPTERS 

In  this  section  we  summarize  the  remaining  chapters  in  this  report. 
In  chapter  2  we  develop  a  general  formulation  for  the  backscattered  echo 
from  a  random  medium  that  can  be  modelled  as  an  ensemble  of  randomly 
distributed  point  scatterers.  To  simplify  the  mathematics  we  assume 
isotropic  scattering.  The  average  backscattered  power  from  the  random 
ensemble  is  calculated  as  a  function  of  transducer  parameters  and  the 
spatial  density  profile  of  the  scatterers.  It  is  shown  that  the  echo 
from  a  random  ensemble  exhibiting  a  sharp  volume  density  gradient  may 
resemble  specular  reflection.  This  explains  for  example  the  "specular" 
reflection  from  organ  boundaries  in  medical  Ultrasound  B  scans.  It  is 
also  shown  that  the  ratio  between  the  power  of  the  "specular"  echo,  l.e. 
the  spatially  coherent  portion  of  the  echo  amplitude  to  the  echo  variance 
should  enable  us  evaluate  the  scattering  density  if  there  is  either  a 
density  gradient  or  a  field  gradient.  The  results  for  the  former  case 
are  verified  qualitatively  by  experiment. 

After  studying  in  chapter  2  the  general  behavior  of  the  backscattered 
echo  from  an  ensemble  of  point  scatterers  we  investigate  in  chapter  3 
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theoretically  an  algorithm  for  clutter  reduction  that  was  first  suggested 
In  [3,4].  The  technique  Is  called  the  "Minimization  Algorithm"  and  In 
fact  provided  the  Impetus  for  this  research  project.  We  calculate  the 
Improvement  In  signal  to  noise  ratio  of  this  algorithm  and  compare  It 
with  experimental  results  obtained  In  [3,4].  The  theoretical  results 
seem  to  fit  well  with  experiment.  Finally  we  calculate  the  Receiver 
Operating  Characteristic  of  this  algorithm. 

In  chapter  4  we  pose  the  problem  of  clutter  reduction  under  the 
assumption  that  the  local  properties  of  the  clutter  echo  are  unknown, 
namely  the  local  variance.  We  describe  a  technique  to  estimate  the  local 
properties  of  the  signal  without  a  priori  knowledge  of  whether  a  target 
exists  In  the  region  of  Interest.  We  use  the  results  to  develop  so-called 
optimum  detection  algorithms  for  both  additive  and  multiplicative  noise 
based  on  some  known  results  In  statistical  decision  theory.  For  the 
Optimum  Detector  algorithm  for  additive  noise,  we  calculate  the  prob¬ 
ability  of  detection  as  a  function  of  signal  to  noise  ratio  (for  the 
clutter  reduction  problem) .  Some  experimental  results  for  this  type  of 
processing  for  the  enhancement  of  flaw  to  grain  echoes  In  metal  are  also 
given.  Resolution  performance  and  range  bias  of  the  suggested  technique 
are  evaluated  experimentally.  We  also  compare  the  performance  of  the 
"Minimization  Algorithm"  with  the  new  technique  on  the  basis  of  ROC 
(Receiver  Operating  Characteristics),  resolution,  bias,  and  amplitude 
dynamic  range  of  the  processed  data. 

Finally  In  chapter  5  we  present  conclusions,  remarks  and  some  sugges¬ 


tions  for  future  work. 


CHAPTER  2 


SCATTERING  FROM  RANDOM  MEDIA 


INTRODUCTION 

In  this  chapter  we  investigate  the  properties  of  the  backscattered 
echo  from  random  media.  The  results  obtained  here  suggest  new  techniques 
for  scatterer  density  estimation  and  also  form  the  basis  for  the  clutter 
reduction  algorithms  to  be  discussed  in  chapter  3  and  4. 

The  random  media  throughout  this  work  are  assumed  to  be  composed  of 
a  high  density  (concentration)  of  point  scatterers  that  can  have  random 
size  (strength)  and/or  random  position.  For  simplicity  we  assume  that  the 
average  distance  between  the  scatterers  is  much  larger  than  their  radius 
thus  allowing  the  assumption  that  their  volume  density  is  Poisson  distri¬ 
buted  (if  their  position  is  random). 

It  is  well  known  that  under  the  assumption  of  plane  wave  approxima¬ 
tion  and  constant  scattering  density  the  backscattered  echo  from  random 
media  can  be  considered  Gaussian  with  zero  mean  provided  the  scatterers 
concentration  is  high  enough  (central  limit  theorem).  However,  if  the 
average  scattering  concentration  (or  size)  is  not  uniform  as  a  function 
of  range  delay,  or  the  sonic  field  exhibits  a  gradient  the  problem  becomes 
more  complex.  Siegert  and  Goldstein  (34)  calculated  the  backscattered 
echo  from  a  continuously  Illuminated  time  varying  random  medium,  composed 
of  Identical  size  scatterers  under  the  assumption  of  plane  wave  illumina¬ 
tion.  They  found  out  that  the  backscattered  echo  should  have  a  non  zero 
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mean  (coherent  echo)  which  Is  proportional  to  the  gradient  of  the  scatter¬ 
ing  density.  Some  years  later  Glotov  (36)  analyzed  the  case  of  a  short 
transmitted  signal  illuminating  a  slab  filled  with  discrete  inhomogenei¬ 
ties  for  both  the  plane  wave  and  spherical  wave  approximation  and  again 
showed  that  the  returned  echo  should  have  a  non  zero  mean. 

In  this  chapter  we  investigate  the  properties  of  the  backscattered 
echo  from  random  media  for  pulsed  transmission  by  extending  the  existing 
theories  to  take  into  account  field  fluctuations  and  the  scattering 
density  profile.  We  start  by  analyzing  the  backscattered  echo  from  a 
regular  lattice  with  equal  spacing  and  random  scattering  cross-section, 
later  we  proceed  to  evaluate  the  returned  echo  from  random  media  with 
uniform  scattering  density  and  finally  we  analyze  the  backscattered  echo 
from  random  medium  with  non  uniform  scattering  density  and  arbitrary 
field  structure. 

THE  REGULAR  LATTICE 

Let  the  echo  detected  by  a  sonic  receiver  at  time  t  after  transmis¬ 
sion  of  a  pulse,  due  to  a  single  point  scatterer  at  f  be  written  as 

Eg(t)  -  PiG(fi,t)  (2.1) 

where  G(f,t)  is  defined  as  -he  impulse  response  of  a  system  composed  of  a 
transmitter/receiver  illuminating  a  scatterer  located  at  and  p  is  a 
proportionality  constant  that  depends  on  the  scatterer  scattering  cross 
section.  Suppose  that  we  illuminate  a  regular  lattice  for  which  the  scat¬ 
tering  cross-section  of  the  particles  is  a  random  variable  (see  Figure  1). 
The  backscattered  echo  from  the  lattice  can  be  written  (assuming  the  Born 


approximation) 
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I  - 


Figure  2.1  The  regular  lattice.  (The  scatterers  within  the  lattice 
were  not  drawn  to  retain  clarity). 


I 
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E(t)  -  I  PiG(?i.t) 
i-1 


(2.2) 


where  a  Is  the  number  of  Che  scacterers  In  the  range  cell,  defined  as 
that  region  In  space  from  which  echoes  are  received  at  time  «.  after 
transmission. 

The  average  power  of  the  backscaCtered  echo  becomes 


P(t)  “I  I  ,t)G(?j,t) 

1-1  j-1  ■' 


(2.3) 


Where  -  means  ensemble  average  (and  not  time  average).  In  order  to 
measure  this  average  one  would  move  the  transducer  with  respect  to  the 
lattice  and  average  the  received  power  for  a  given  range  delay  t.  As¬ 
suming  PI  to  be  statistically  Independent  of  pj  eq.  (3)  becomes 


P(t)  -  ^  |G(?i.t)|2  +  p2  ^  I  G(?i,t)G(?j,t) 

1-1  i-1  j-1 


Rewriting  eq.  (4)  we  obtain 


(2.4) 


P(t)  -  I  |G(?i,t)|2  +  p2  {G(?,  ,t)IG(f,  ,t)  + 
1-1 

+  G(f2.t)  +  ...+  G(fn,t)  -  G(?^,t)]  +  +  } 


(2.5) 


Replacing  the  summations  by  Integrals  we  obtain 


fW  -  ii  /  |G(f,t)l2dv  +  IL  1/  G(lf,t)dv|2 

Av  V  ^yZ  V 


-  P_  /  |G(f,t)|2dv 

Av  V 


(2.6) 


where  Av  is  the  volume  of  a  unit  cell  of  the  lattice,  but 


_L_  -  £  -  N 
Aw  V 


(2.7) 
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and*  eq.  (6)  becomes 

P(t)  -  /  lG(f,t)l2dv  +  N2p2  1/  G(f,t)|2 

V  V 

-  NT52  /  lG(<’,t)[2dv  (2.8) 

V 

as  we  see  in  the  next  section  the  second  term  In  eq.  (8)  very  often 
vanishes.  In  this  case  P(t)  becomes 

P(tT  -  Nap2/  |G(?.t)l2dv  (2.9) 

and  the  backscattered  power  is  directly  proportional  to  the  variance  of 
the  scattering  cross-section  of  the  lattice  particles,  o^.  It  Is  obvious 
that  If  the  scatterers  are  uniform  In  size  no  backscattered  power  Is 
expected. 

RANDOM  MEDIUM  WITH  UNIFORM  SCATTERING  CONCENTRATION 

When  a  random  ensemble  Is  Illuminated  by  a  transducer  the  received 
echo  at  time  t  after  transmission  can  be  written** 

n 

E(t)n)  -  I  PiG(?i,t)  (2.10) 

1-1 

where  n  Is  the  number  of  the  scatterers  In  the  range  cell.  We  assume  that 

the  range  cell  has  a  volume  V,  and  that  It  contains  n  randomly  positioned 

♦  ♦  +  -*■ 

scatterers  at  locations  ri,r2,r3, . . . .r^  of  strengths  Pi, P2 > • • • » Pn* 
Instantaneous  power  at  time  t  Is 

*  V 

For  non  cubic  unit  cell  ^wlll  be  directly  proportional  to  n  and  not 
equal  to  n  as  Is  assumed  here  for  a  cubic  unit  cell. 

**E(t|n)  should  be  read;  E  at  range  delay  t  given  that  there  are  exactly 
n  scatterers  In  the  range  cell. 
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a  n 

P(t|Q)  -  I  I  Pi  Pi  G(?.,t)  (2.11) 

1-1  j-1 

This  can  be  separated  Into  two  parts 

P(t|n)  -  I  |G(?i,t)l2p2  +  5  5  G(fi,t)G(?j,t)PiPj 

1-1  i  1-1  j-1  (2.12) 

i^j 


The  average  power  at  time  t  after  transmission,  for  differing  locations 
of  the  n  scatterers  In  the  range  cell  of  volume  V  can  be  written 


P(t|n)  -  n  p^|G(?jL»t)P  +  n(n-l)  pj^  p^  G(?j^,t)  G(f^,t) 

■'  ■'  (2.13) 

Note  that  P(t|n)  Is  the  average  over  all  possible  configurations  of  the 

scatterer  location  In  the  range  cell.  Assuming  that  Is  uncorrelated 

with  f'j  for  l^j,  we  can  rewrite  eq.  (2.13)  as 

P(t|n)  -  n  /  lG(^,t)l2  p(^)dv  +  n(n-l)  p^  |  /  G(f,t)p(f’)dvl2 

V  V  (2.14) 

where  p(f’)  Is  the  probability  of  finding  one  Individual  scatterer  In  the 
voliune  element  at  t,  and  where  the  Integrals  are  taken  over  the  volume  V 
of  the  range  cell.  We  may  write 

P(if’)  -  -  (2.15) 

V 

Using  these  two  relations  we  can  write  the  ensemble  averaged  power  P( t | n) 

P(tln)  -  -  p2  /  lG(lf,t)l2dv  +  n(nilip2  |  /  G(lf’,t)dv|2  (2.16) 

V  V  V2  V 


To  find  the  average  power  P(t)  we  have  to  evaluate 
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P(t)  -  I  P(tin)  p(Q)  (2.17) 

n-0 

from  eq.  (16)  and  (17) 


P(t)  -  /  |G(f,t)|2dv  I  np(n)  + 

V  V  n-1 


+  1/  G(f’,t)dv|2  I  a(n-l)p(n)  (2.18) 

V2  V  n-1 


where  n  Is  Poisson  distributed  with  n  as  Its  average.  Replacing 


I  np(n)  -  a  (2.19) 

n»0 

and  (as  Is  known  for  the  Poisson  distribution). 


I  n(n-l)  p(n)  -  h2  (2.20) 

n««0 


we  obtain 

p(tT  -  N  p?  /  lG(f,t)l2dv  +  n2  1/  G(f,t)dvl2  (2.21) 

V  V 


where 

a 

N  -  - 
V 


Note  again  that  P(t)  Is  ensemble  average  and  not  time  average  and  that  to 
obtain  this  average  one  would  have  to  move  the  transducer  with  respect  to 
the  medium  (If  the  medium  Is  stationary)  and  average  the  received  power 
for  a  given  t.  He  see  that  the  first  term  on  the  right  hand  side  of  eq. 
(21)  Is  proportional  to  N^  l.e.  to  the  average  of  the  sum  of  the  In¬ 
dividual  backscattered  powers.  We  will  therefore  refer  to  It  as  the 
Incoherent  power  P^nc'  second  term  Is  proportional  to  (Np)  ,  l.e.  to 
the  square  of  the  sum  of  the  backscattered  amplitudes.  We  will  therefore 
describe  It  as  the  coherent  backscattered  power,  Pcoh* 
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We  now  show  that  P^oh  equals  the  magnitude  squared  of  the  averaged 
echo,  and  that  Pine  equals  the  variance  of  this  echo. 

From  eq.  (10),  (15)  and  (19) 


E(t)  -  I  E(t|n)  p(n)  -  Np  /  G(f, t)dv 
n-0  V 


(2.22) 


So  that 


Ie|2  -  n2  f  1/  G(?,t)dv|2 

V 


(2.23) 


The  fact  that  Pine  the  variance  of  E(t)  follows  Immediately,  since 


4  -  lEp  -  1e|2 


(2.24) 


Thus  from  eqs.  (21)  and  (24), 


®E  “  ^inc  -  /  lG(?,t)l2dv 

V 


(2.25) 


From  eq.  (21)  one  can  see  that  the  coherent  term  will  be  appreciable  if 
either  N  or  /  G(f,t)dv  are  sufficiently  large.  Let  us  examine  the  Inte- 

V 

gral  /  G(f,t)dv  for  simple  beam  geometries.  For  example  consider  a  plane 

V 

wave  travelling  along  the  z  axis  for  which  the  field  function  can  be 
written 


G(f,t)  -  I(t  -  2z)dz 
c 


(2.26) 


The  coherent  term  of  eq.  (22)  becomes 


E^  -  Np  /  I(t  -  2z)dz 

.flB  C 


(2.27) 


Changing  variables  t'  -  t-2z/c  eq.  (22)  becomes 


irry  -  /  I(t')dt’ 

2  rnrnm 


(2.28) 
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The  Integral  and  with  It  the  coherent  term  reduces  to  zero  since  the 
signal  cannot  contain  a  dc  component.  Likewise  for  a  point  transmitter/ 
receiver 

Oft 

E(tT  -  Np  /  (I(t-2r/c)/r2)4Trr2dr  (2.29) 

0 

which  again  reduces  to  zero.  Thus  we  see  that  both  for  the  case  of  plane 
wave  approximation  as  well  as  for  the  case  of  waves  emitted  from  a  point 
transmitter/receiver  the  coherent  term  vanishes. 

However  in  a  work  which  is  currently  performed  by  Goyao  Yu  from 
Drexel  University  it  is  shown  that  for  almost  any  other  geometry  of  trans¬ 
mitter/receiver  the  integral  /  G(f, t)dv  does  not  vanish.  For  example, 

V 

Mr.  Yu  calculated  numerically  this  Integral  for  a  point  transmitter  and 
a  ring  shaped  receiver  and  showed  that  the  coherent  integral  is  different 
from  zero  in  the  near  field.  However  In  the  far  field  the  coherent  inte¬ 
gral  practically  vanishes.  It  is  clear  that  the  Integral  |/  G(f,t)dv|2  is 

V 

usually  much  smaller  than  /  |G(f,t)|2dv.  However  the  Integral  |/G(f,t)dv|2 

V  V 

is  multiplied  by  while  the  integral  /  |G(f,  t)|2<iv  is  multiplied 

V 

only  by  N  (in  eq.  2.21)  which  should  make  the  ratio  |El^/a2g  measurable. 

These  results  obtained  by  Mr.  Yu  give  rise  to  the  hope  that  scatterer 
concentration  and  scattering  cross  section  estimation  would  be  possible 
using  coherent  reflection  caused  by  field  gradients. 

RANDOM  MEDIUM  WITH  NONUNIFORM  SCATTERING  DENSITY 

If  the  scatterers  exhibit  a  certain  scattering  concentration  profile 
N(f)  the  algebra  introduced  earlier  becomes  somewhat  more  Involved.  Rede¬ 
fining  pCi?’)  to  be  N(if’)i.  the  average  power  at  a  certain  range  delay  t  can 
be  shown  to  be, 
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P(t) 


/  N(f)lG(?,t)l2dv  + 
V 


_2 

P 


1/  N(f)  G(r,t)dvl2 
V 


(2.30) 


with 

ECO  “  P  /  N(f)G(t,t)dv  (2.31) 

V 

It  can  be  seen  by  Inspection  that  the  Integral  of  eq.  (31)  will  be  non¬ 
zero  provided  N(0  has  a  gradient  In  the  same  direction  as  that  of  the 
function  G(f,t),  l.e.  In  the  direction  of  the  sound  beam.  For  the  Inte¬ 
gral  to  be  large  the  gradient  of  N(f)  should  be  sharp  compared  to  the 
sound  wavelength.  The  value  of  the  Integral  for  a  step  function  In  density 
N(0  Is  derived  below.  Eq.  (31)  shows  that  a  gradient  In  scattering 
density  In  the  direction  of  propagation  of  the  sound  beam  will  return  a 
spatially  coherent  echo.  The  same  holds  true  for  the  echo  from  the  bound¬ 
ary  between  two  regions  with  different  scattering  densities.  For  simplic¬ 
ity  we  assume  the  Imaginary  boundary  to  be  planar.  Replacing  f  by  z  we 
can  write, 


N(z)  -  z  <  Zq  -  (2.32) 

-  N2  z  >  Zq 

For  a  plane  wave  G(z,t)  ■  b(t-2z/c)cos[u(t-2z/c)  +  9];  thus  from  eq.  (31) 
the  returned  echo  from  the  boundary  becomes 

m  O 

E(tQ)  -  £  pj  /  b(t)co8[  a)t+9]dt  +  £  N2  P2  /  b(t)co8  [  u>t+9] dt 
2  Q  2  —CD 

(2.33) 

and  P2  3re  the  average  reflection  coefficients  of  the  scatterers  In 
region  1  and  2  respectively.  Note  that  the  volume  Integral  In  eq.  (30) 
becomes  a  time  Integral  by  change  of  variables .  If  the  range  cell  Is 
situated  such  that  the  boundary  Is  at  Its  center  the  variance  of  E(tQ)  can 


be  bliown  t J  be , 
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a2  H  i[o?  +  a^]  3  £  [N,"^  +  No"^]  /  lb(t)l2dt 

2  ^  4  ^  ^  o  (2.34) 

2 

where  a|  and  a2  are  the  variances  of  the  returned  echo  from  region  1  and 
region  2  respectively.  For  N2“0  eqs.  (33)  and  (34)  become 

OD 

E(tQ)  -  Nj  -  /  b(t)co8[a)t+0]dt  (2.33a) 

2  o 


SNELL'S  law  for  DENSITY  GRADIENTS 

It  Is  easy  to  show  that  the  scattered  sound  from  density  gradients  In 
a  random  scatterer  ensemble  obeys  Snell's  law.  Consider  two  plane  wave 
transducers  used  for  transmission  and  reception  respectively,  oriented 
at  +  9  and  -  9  degrees  with  repect  to  the  z  axis.  Then  the  <  and  f 
components  of  7  G(f,t)  are  zero  where  the  function  G(f,t)  Is  as  defined 
In  eq.  (1)  but  for  separate  transmitter  and  receiver. 

Although  eq.  (31)  above  was  derlveo  for  the  case  of  a  single  trans¬ 
ducer  used  for  both  transmission  and  reception,  It  Is  equally  valid  for 


the  transducer  pair  described  In  the  previous  paragraph.  For  plane  wave 
transmission  and  reception  It  can  be  shown  that  the  Integral  In  eq.  (31) 
Is  non  zero  provided  that  the  gradients  of  N(l^)  and  G(l^, t)  are  parallel. 
Thus  a  density  gradient  In  a  random  scatterer  ensemble  will  scatter  a 
coherent  Incident  beam  In  the  direction  given  by  Snell's  law,  l.e.  In  the 
same  direction  In  which  the  beam  would  be  reflected  by  a  mirror  parallel 
to  the  planar  constant  density.  This  result  Is  not  unexpected  and  Is 
encountered  often  (without  much  excitement)  In  almost  every  field  of 
ultrasonic  measurement. 

EXPERIMENTAL  RESULTS 

According  to  eq.  (33)  we  should  be  able  to  observe  a  coherent  echo 
from  the  boundary  between  two  different  regions  of  differing  scattering 
density. 

One  of  the  easiest  ways  to  obtain  a  sharp  boundary  for  scatterers 
Immersed  In  water,  Is  to  use  sponges.  A  sponge  can  be  cut  to  produce  a 
sharp  boundary  and  when  Immersed  In  water  the  complex  structure  of  the 
sponge  fibers  can  be  considered  as  randomly  distributed  scatterers. 
Furthermore,  one  can  clamp  together  two  sponges  with  different  scattering 
properties  and  thus  compose  a  sharp  boundary  between  two  different  media, 
or  one  can  use  only  one  sponge,  simulating  two  media  with  N^-O  In  region 
1  and  N2  In  region  2.  Figures  2(A)-2(D)  show  pictures  of  sponges  used  In 
such  experiments.  One  can  see  that  sponge  A  has  the  finest  honeycomb 
structure.  Sponge  B  has  a  larger  honeycomb  structure,  sponge  C  Is  even 
more  dilute  and  sponge  D  is  the  most  dilute  In  comparison  to  sponges  A,B 
and  C.  In  Figure  3  we  see  sponges  B  and  C  side  by  side.  One  can  see 
that  the  boundary  Is  sharp  In  comparison  to  a  wavelength  (frequency  of 


Figure  2.2.  A  photograph  of  the  surfaces  of  sponges  A,B,C  and  D. 
fine  honeycomb  structure  Is  visible.  (The  width  of  each  picture 
responds  to  6.5  mm.) 


Figure  2«3«  Sponges  B  end  C  side  by  side*  It  Is  seen  thet  the  honeycomb 
structure  of  sponge  B  Is  finer  than  that  of  sponge  C> 


^  ^  \A  At- 


Figure  2.4.  Four  different  echoes  from  sponge  B.  Note  that  the  coherent 
echo  component  from  the  sponge  boundary  and  the  Incoherent  component  from 
within  the  sponge. 
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2.25  MHz).  One  can  also  see  again  that  the  honeycomb  structure  of  sponge 
B  is  finer  than  that  of  sponge  C. 

Figure  4  shows  echoes  returned  from  sponge  B.  One  can  see  that  the 
echoes  from  the  boundary  are  highly  coherent  and  exhibit  negligible  phase 
change  In  comparison  to  the  more  random  echoes  from  Inside  the  sponge. 
These  echoes  are  Incoherent  in  the  sense  that  for  a  certain  range  delay 
the  amplitude  la  random  with  zero  mean.  Note  that  the  echoes  In  close 
vicinity  to  the  first  echo  seem  somewhat  weaker  than  those  from  deeper  in¬ 
side  the  sponge  (at  greater  depths  the  echoes  weaken  due  to  attenuation). 
The  reason  Is  probably  due  to  the  fact  that  the  transducer  surface  Is 
not  perfectly  parallel  to  the  sponge  surface  and  thus  It  takes  more  time 
for  the  Incoherent  term  to  fully  develop.  This  slight  misorientation  has 
practically  no  effect  on  the  coherent  term  magnitude.  However  because  of 
this  effect,  the  estimation  of  the  Incoherent  term  should  not  be  done  In 
the  Immediate  vicinity  of  the  first  echo. 

From  eq.  (35)  we  know  that  the  ratio  between  the  square  of  the  aver¬ 
age  boundary  echo  and  the  average  power  from  Inside  the  sponge,  neglecting 
attenuation.  Is  proportional  to  the  scattering  density.  Thus  It  should 
be  possible  to  estimate  the  spouge  density  from  A-mode  echoes  of  the 
types  obtained.  We  will  not  attempt  here  to  estimate  the  sponge  density 
quantitatively,  but  will  show  qualitatively  that  the  experimental  results 
behave  according  to  the  theory. 

Figures  5(i'  .(B), (C)  and  (D)  represent  the  echoes  received  from 
sponges  A,B,C  and  D  respectively.  Each  of  these  pictures  Is  composed  of 
four  superimposed  traces  from  four  different  locations  In  the  sponge, 
with  Identical  distances  between  the  transducer  and  the  sponge  surface. 
In  Figure  5(A)  which  represents  the  received  echo  from  sponge  A  one  can 
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: 

<  -- 

A 

Figure  2.5(A)  Four  traces  (one  on  top  of  the  other)  of  the  reflected  echo 

from  sponge  A.  Note  the  large  coherent  component  of  the  echo  from  the 

boundary. 

(B)  Sponge  B.  Note  that  the  coherent  component  is  not  as  large  in 

comparison  to  the  echo  from  within  the  sponge  as  in  Figure  4(A) 

(C)  Sponge  C.  Note  that  the  coherent  component  is  smaller  in  comparison 
to  the  echo  from  within  the  sponge  in  Figure  4(B) 

(D)  Sponge  D.  The  coherent  component  was  found  to  be  negligible  in 

repetitive  experiments. 
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see  clearly  that  the  coherent  echo  is  large  In  comparison  to  the  echo 
returned  from  within  the  sponge.  Figure  5(A)  is  similar  to  Figure  4  but 
with  the  traces  superimposed.  One  can  see  that  the  coherent  component 
of  the  echo  from  the  boundary  is  not  as  large  as  that  of  Figure  5(A)  in 
comparison  to  the  echo  returned  from  within  the  sponge. 

Notice  also  the  great  difference  in  the  degree  of  coherency  for 
sponges  B  and  C.  The  honeycomb  structures  of  these  sponges  are  similar  in 
shape  so  that  we  can  assume  tha':  ^2/^  is  also  similar.  The  integrals 
in  eq.  (30)  are  also  similar  for  the  two  sponges  since  they  involve  only 
transducer  and  medium  parameters  which  were  the  same  for  all  experiments. 
Thus  the  ratio  between  the  coherent  term  and  the  incoherent  term  of  sponge 
B  and  C  should  be  proportional  to  the  scatterer  density.  Evaluating  the 
coherent  term  from  the  ratio  of  the  average  of  the  first  echo  from  the 
boundary  to  the  power  from  inside’  the  sponge,  (the  incoherent  term)  we 
find  the  ratio  to  be  about  2.1  for  sponge  B  and  about  0.8  for  sponge  C 
which  implies  a  density  ratio  of  about  2.6.  The  actual  scatterer  density 
that  can  roughly  be  established  from  the  micrographs  is  about  2.2  which 
is  in  sufficiently  close  agreement  considering  the  fact  that  only  four 
sample  points  were  used  and  that  no  special  arrangements  were  made  for 
producing  extremely  smooth  surfaces  and  for  keeping  the  sponge  surface 
parallel  to  the  transducer  surface.  Figure  5(D)  corresponding  to  sponge  D 
with  the  largest  honeycomb  structure  has  practically  zero  coherent  effect. 
This  is  due  to  the  fact  that  sponge  D  has  much  lower  density  than  the 
other  sponges,  too  low  to  show  a  coherent  term  from  only  4  sample  points. 

Figure  6  shows  the  echo  from  sponge  D  clamped  to  sponge  B.  Coherent 
echoes  are  clearly  seen  along  the  center  vertical  line  of  the  picture. 

The  angle  dependence  of  the  reflected  echo  was  also  investigated, 


Figure  2.6.  The  received  echo  from  a  2  layer  sponge  complex  composed  of 
sponges  D  and  B.  The  coherent  component  can  be  seen  on  the  central 
vertical  line-  Also  note  that  the  power  reflected  from  sponge  D  is  less 
than  that  reflected  from  sponge  B. 


Figure  2.7.  Angle  dependence  of  the  coherent  component  (sponge  B). 
Note  that  while  the  coherent  term  is  highly  sensitive  to  angle  variations, 
the  echo  from  within  the  sponge  is  angle  insensitive. 
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using  sponge  B  as  a  reflector.  One  can  see  from  Figure  7  that  the  coherent 
effect  Is  highly  angle  sensitive  whereas  the  echo  reflected  from  within 
the  sponge  Is  not  angle  dependent.  In  fact  the  behavior  of  the  coherent 
component  Is  very  much  like  that  of  a  specular  reflection  as  predicted  by 
the  theory.  The  reflection  from  within  the  sponge  Is  Independent  of  angle 
since  this  echo  must  be  Independent  of  boundary  region. 

It  Is  pertinent  to  Inquire  how  well  the  sponges  used  as  test  objects 
In  this  work  fulfill  the  requirements  of  the  theoretical  model.  This 
model  places  no  restriction  on  the  size  of  the  Individual  scatterers  or 
on  the  number  of  the  scatterers  In  the  rauge  cell  but  assumes  the  Poisson 
distribution  which  Implies  that  the  scatterers  are  separated  by  a  distance 
large  compared  to  their  diameter.  Weak  scattering  Is  also  assumed  which 
Implies  that  each  scatterer  produces  no  more  than  one  echo.  Figure  2  in¬ 
dicates  that  the  size  of  the  sponge  unit  cells  Is  no  larger  than  that  of 
the  sound  wavelength  of  about  0.6  mm.,  but  does  not  provide  Information 
about  the  size  separation  of  the  Individual  scatterers.  The  fact  that 
the  scatterers  echoes  from  within  the  sponge  are  found  to  be  spatially 
Incoherent,  shows  that  the  scatterer  distribution  even  In  the  hexagonal 
lattice  sponges  B  and  C,  Is  at  least  pseudo-random.  To  establish  the  fact 
that  the  coherent  reflection  Is  not  due  to  a  change  In  acoustic  Impedance, 
we  measured  the  sound  velocity  within  the  Immersed  sponges,  and  found  It 
to  be  similar  to  the  sound  velocity  In  water. 

DISCUSSION 

As  can  be  seen  In  each  of  the  analyzed  cases  the  ratio  between  the 
coherent  term  and  the  Incoherent  term  for  the  plane  wave  approximation 
Is  proportional  to  the  scattering  density  profile  and  the  function  G(r,t) 
(eq.  (35)).  In  order  that  the  upper  Integral  will  be  appreciable  the 


31 


scattering  density  gradient  should  be  sharp  compared  to  the  wavelength. 
If  the  gradient  is  not  sharp  enough  the  Integral  and  with  it  the  coherent 
term  become  negligible  and  the  backscattered  echo  can  be  considered  as 
Gaussian  (if  the  scattering  density  is  high  enough)  with  zero  mean  and 
variance  which  is  a  function  of  the  range  delay.  However,  if  the  gradient 
is  sufficiently  large  the  coherent  term  may  be  appreciable  even  if  the 
scattering  density  is  low. 

It  was  shown  that  the  reflection  from  a  sharp  boundary  separating 
regions  of  different  scattering  density  (or  scattering  cross-section)  is 
characterized  by  having  a  non  zero  mean  when  spatially  averaged  and  a 
variance  of  approximately 

i  [d^  +  a|]  (2.34) 

The  statistics  of  the  boundary  echo  can  be  considered  Gaussian  due 
again  to  the  central  limit  theorem  with  mean  E(t)  (eq.  (33))  and  variance 
according  to  eq.  (34).  This  result  states  in  fact  that  if  the  scattering 
density  gradient  is  sharp  enough  the  reflection  from  the  boundary  between 
regions  of  different  scattering  densities  resembles  specular  reflection 
and  it  might  be  appropriate  to  model  the  boundary  as  a  strong  specular 
reflector  embedded  in  a  high  density  of  point  scatterers.  We  can  use 
this  result  to  Incorporate  boundary  detection  in  chapters  3  and  4  where  we 
analyze  some  detection  schemes  for  targets  embedded  in  clutter  environ¬ 
ments  . 

We  also  showed  theoretically  that  the  ratio  of  the  coherent  signal 
from  the  boundary  of  a  scattering  region  to  the  incoherent  signal  from 
Inside  the  region  (if  the  scattering  density  is  low  enough)  is  propor¬ 
tional  to  the  average  spatial  density  of  the  scatterers  multiplied  by  a 
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shape  factor.  The  experimental  results  suggest  that  this  technique  might 
be  useful  In  nondestructlvely  estimating  scatterer  size  and  density  even 
If  the  scatterers  are  much  too  small  and  much  too  close  to  be  resolved. 
It  might  also  be  applied  to  extract  more  Information  on  the  tissue  state 
In  biomedical  ultrasonic  Imaging  In  regions  for  which  the  boundary  struc¬ 
ture  Is  known  and  Is  smooth  compared  to  the  Illuminating  wavelength. 

Another  practical  application  that  this  method  might  provide  Is  In 
the  estimation  of  the  surface  grain  size  of  metals.  In  order  for  this 
technique  to  succeed  the  metal  would  have  to  be  polished  and  Immersed  In 
a  liquid  of  similar  characteristic  Impedance. 


CHAPTER  3 


I 


THE  MINIMIZATION  ALGORITHM 


INTRODUCTION 

In  chapter  1  we  outlined  several  algorithms  and  techniques  for  the 
reduction  of  clutter  and  speckle.  We  found  out  that  the  techniques  In 
most  cases  do  not  provide  sufficient  Improvement  In  clutter  reduction. 
The  need  for  post  reception  algorithms  for  the  enhancement  of  target  to 
clutter  ratio  Is  evident.  In  chapters  3  and  4  we  will  Investigate  two 
algorithms  for  the  reduction  of  clutter  using  split  spectrum  processing. 
This  type  of  processing  and  the  algorithm  to  be  described  In  chapter  3 
was  first  Introduced  by  Newhouse  et  al  in  1982  (3).  N.  Bllgutay  in  his 
Fh.D.  dissertation  empirically  Investigated  this  algorithm  the  so-called 
"Minimization  Algorithm",  and  showed  experimentally  that  it  provides  an 
Improvement  In  slgnal-to-nolse  ratio  over  conventional  techniques.  In 
this  chapter  we  analyze  this  algorithm  theoretically.  We  evaluate  its 
slgnal-to-nolse  ratio  and  its  Receiver  Operating  Characteristics  and 
compare  Its  performance  to  other  split-spectrum  algorithms. 

Chapter  4  Introduces  a  new  algorithm  which  we  call  "Optimal  Pro¬ 
cessing"  (The  processing  Is  optimal  In  the  sense  that  for  every  range 
delay  the  processor  maximizes  the  detection  performance). 

These  two  chapters  can  be  read  independently.  Historically  the  Mini¬ 
mization  Algorithm  proceeded  the  Optimal  Processing  and  was  invented  "by 
Intuition".  So  we  start  with  the  analysis  of  the  Minimization  Algorithm 


and  procceed  In  chapter  4  to  Introduce  the  Optimal  Processor.  We  compare 
the  two  algorithms  on  the  basis  of  detection  properMes  at  the  end  of 
chapter  4. 

The  split  spectrum  processing  was  Introduced  by  V.  L.  Newhouse,  N.  M. 
Bllgutay  and  E.  S.  Furgason  [3]  In  1979.  In  this  technique  the  wideband 
spectrum  of  a  received  signal  is  split  into  several  frequency  windows. 
The  resulting  signals  are  Independent  provided  that  these  frequency 
windows  do  not  overlap,  and  can  be  compounded  In  ways  similar  to  those 
used  in  frequency  agility  so  as  to  reduce  clutter  noise  relative  to  target 
echo  (see  Figure  1).  Following  the  nomenclature  Introduced  by  Abbot  and 
Thurston  [19]  for  their  technique  of  Illumination  with  different  frequen¬ 
cies,  we  will  refer  to  their  technique  and  to  that  of  reference  [3]  as 
"frequency  compounding"  so  as  to  bring  out  the  similarity  between  these 
techniques  and  that  of  spatial  compounding. 

Of  the  frequency  compounding  techniques  described  by  Newhouse  et  al 
[3],  the  most  powerful  used  a  new  procedure,  (the  so-called  Minimization 
algorithm).  In  which  the  signals  obtained  by  frequency  splitting  the  echo 
spectrum  are  squared  and  the  minimum  value  at  each  delay  Is  chosen.  This 
technique  was  found  to  provide  a  significant  Improvement  In  slgnal-to- 
nolse  ratio  with  respect  to  the  compounding  techniques  previously  used 
with  frequency  agility  In  either  radar  or  ultrasound  in  which  the  signals 
obtained  by  frequency  splitting  were  either  averaged  and  then  squared  or 
were  first  squared  and  then  averaged.  A  similar  technique  has  recently 
been  applied  to  two-dimensional  imaging  [41]. 

This  chapter  Is  devoted  to  a  theoretical  analysis  of  the  Minimization 
Algorithm  and  Its  comparison  with  the  two  earlier  known  frequency  agility 
algorithms.  We  present  expressions  for  signal  to  noise  ratio,  probability 
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of  detection  and  false  alarm  as  a  function  of  the  number  of  frequency 
windows,  target  size  and  ratio  between  the  clutter  variance  within 
and  outside  the  range  cell.  The  analytical  results  are  compared  with 
simulated  zvA  expsrlmen*-al  data  tak?n  from  (4).  Finally  it  is  shown 
analytically  that,  In  most  cases,  the  Minimization  algorithm  provides 
a  much  larger  slgnal-to-nolse  ratio  enhancement  than  the  two  classical 
algorithms.  However,  It  Is  found  that  this  Improvement  In  slgnal-to-nolse 
ratio  enhancement  results  In  a  trade-off  which  Increases  the  probability 
of  error  in  detection. 

Although  the  algorithms  described  here  could  have  applications  In 
several  different  fields.  It  should  be  noted,  that  the  assumptions  made 
In  the  analysis  refer  only  to  ultrasonic  Imaging,  In  either  medical 
applications  or  nondestructive  evaluation. 


ANALYSIS 

In  this  section  we  perform  a  comparative  analysis  of  the  three  fre¬ 
quency  compounding  algorithms  under  consideration.  We  first  define  our 
model  for  the  target  and  Its  surrounding  randomly  distributed  scatterers, 
and  specify  our  definition  of  signal  to  noise  ratio.  We  then  compute  the 
signal  to  noise  ratio  enhancement  of  the  three  algorithms  under  considera¬ 
tion,  and  finally  calculate  their  Receiver  Operating  Characteristics  (ROC). 

In  all  three  of  the  algorithms  analyzed  here  we  obtain  the  frequency 
diverse  signal  set  used  for  the  frequency  compounding,  by  filtering  (or 
"splitting”)  the  received  wide  band  echo  from  the  target-clutter  source 
complex  Into  different  frequency  bands  or  windows  which  may  or  may  not 
overlap.  For  the  sake  of  mathematical  simplicity  our  analysis  In  this 
chapter  Is  limited  to  the  case  of  non-overlapping  frequency  windows. 


Note  that  the  length  of  the  range  cell  of  the  frequency  compounded  image 
Is  Inversely  proportional  to  the  bandwidth  of  these  frequency  windows. 

Target  Model  and  Signal  to  Noise  Ratio 

The  model  to  be  analyzed  Is  composed  of  a  strong  point  reflector  (the 
target)  embedded  In  a  high  density  of  smaller  scatterers  (the  clutter). 
We  recognize  two  differences  between  the  properties  of  the  target  and  any 
of  the  scatterers. 

1.  The  backscattered  echo  from  the  target  Is  much  larger  than  that  of  a 
point  scatterer. 

2.  There  is  only  one  target  In  the  range  cell. 

In  order  to  determine  the  statistics  of  the  peak  target  echo  in  the 
presence  of  the  clutter  we  perform  a  Gedanken  experiment  In  which  we  keep 
the  target  In  Its  place  but  change  the  clutter  configuration  from  experi¬ 
ment  to  experiment  and  measure  the  echo  voltage  at  a  range  delay  corre¬ 
sponding  to  the  peak  target  echo.  The  statistics  of  the  clutter  without 
a  target  In  the  range  cell  will  be  Gaussian  (see  for  example  (1)).  It 
follows  that  the  statistics  of  the  echo  corresponding  to  the  target  range 
delay  will  also  be  Gaussian  but  with  non-zero  mean.  Thus  we  could  write 
two  hypotheses  for  this  case 

Hi;  r  *  m  +  ni  for  the  case  of  the  range  cell  containing  noise  (3.1) 
and  a  target 

Hq:  r  >  ni  for  the  case  of  the  range  cell  containing  noise  only 

Where  r  Is  the  value  of  the  received  signal  at  the  analyzed  delay  t,  m  Is 
defined  as  the  target  peak  echo  and  n^  Is  the  noise  voltage  which  Is  a 
superposition  of  the  contribution  of  the  scatterers  (clutter)  In  the  range 
cell,  n^  Is  a  zero  mean  Gaussian  clutter  noise  with  standard  deviation 
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(s.d.)  If  we  knew  the  statistical  parameters  of  the  signal,  (i.e. 

and  m)  the  statistical  decision  rule  would  be  trivial.  For  example  set¬ 
ting  a  threshold  level  at  m/2  would  provide  the  best  decision  rule  for 
minimum  probabillt./  of  error.  In  ultrasonic  imaging  however  we  do  not 
have  a  priori  information  on  either  m  or  a][.  Furthermore  these  parameters 
change  as  a  function  of  distance  from  the  transducer  due  to  attenuation, 
shadowing  or  fluctuations  in  scattering  concentration  or  scattering  cross 
section.  In  order  to  evaluate  the  different  techniques  (to  be  presented 
shortly)  we  introduce  the  hypothesis  set  in  the  following  ways 

r  *  m+ni 
Hq:  r  -  no 

Here  n^  and  ng  are  zero  mean  Gaussian  clutter  noise  echoes  with  s.d  a]_ 
and  Oq  respectively;  ng  will  be  considered  always  to  have  equal  or  larger 
variance  than  n^  (m  is  as  before).  The  need  for  introducl  different 
noise  levels  outside  the  range-cell  is  to  enable  correct  evaluation  of 
the  minimization  algorithm.  This  algorithm  is  non-linear  and  thus  treats 
signal  plus  noise  differently  from  noise  alone  (depending  on  the  local 
slgnal-to-noise  ratio).  Therefore  with  this  model  we  are  able  to  test  the 
effect  of  algorithms  on  clutter  with  unknown  parameters  that  may  slowly 
vary  as  a  function  of  position.  To  complete  the  set  we  should  have  added 
a  third  hypothesis,  H2 

H2:  r  -  03 


for  the  case  of  noise  far  away  from  the  target  which  is  smaller  than  that 
near  the  target.  For  simplicity  we  exclude  this  hypothesis  from  the  set. 
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This  will  be  shown  to  be  Justified,  since  Figure  3  demonstrates  that  the 
algorithms  considered  here  are  monotonlc  In  the  sense  that  If  n3<nQ  then 
the  output  of  the  processed  noise  03  Is  smaller  than  the  processed  noise 
Hg.  Thus  enhancement  of  over  Hg  la  a  sufficient  measure  of  the  Improve¬ 
ment  In  the  performance  of  the  algorithm. 

Since  the  signal  must  be  squared  for  two  of  the  three  algorlt.ims  con¬ 
sidered  below,  we  will  define  the  Input  slgnal-to-nolse  ratio  In  terms 
of  the  square  of  the  Input  signal  (unspllt  broad  band  echo) .  Even  this 
condition,  however.  Is  not  sufficient  to  define  SNR  unambiguously.  We 
therefore  select  a  definition  which  matches  as  closely  as  possible  the 
manner  In  which  SNR  Is  computed  In  the  simulations  and  experiments  with 
which  our  analysis  Is  compared  below. 

We  denote  as  the  square  of  the  Input  signal  for  case  (target 
In  the  range  cell),  and  n^n  as  the  square  of  clutter  signal  amplitude  for 
case  Hg.  We  define  the  SNR  as 


^  hln2 


(3.2) 


The  average  of  the  Input  signal  defined  as  the  square  of  the-  target  echo 
Is: 


Zln  “  -  E((m  +  ni)^}  -  m^  +  <,^2  (3.3) 

and  the  clutter  noise  outside  the  target  region  Is 

-  [E(no^}  .  [30g^]^/2  .  /3-  ^^2  (3.4) 

Thus  the  input  SNR  becomes: 
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2  2 
m  +02^ 

/3  a2 


(3.5) 


Signal  to  Noise  Ratio  Enhancement  for  Non-Linear  Averaging 

We  are  now  in  a  position  to  calculate  the  signal-to-noise  ratio  en¬ 
hancement  of  the  so-called  non-linear  averaging  algorithm,  by  which,  when 
used  In  radar  frequency  agility  systems.  Incoming  signals  at  different 
frequencies  are  first  squared  and  then  averaged.  When  used  In  frequency 
compounding,  as  Is  to  be  analyzed  here,  we  start  by  splitting  the  received 
wide  band  echo  Into  n  adjacent  but  not  overlapping  frequency  bands  or 
windows.^  This  procedure  Is  equivalent  to  transmitting  the  frequency 
bands  consecutively  as  is  done  In  frequency  agility.  Therefore  the  results 
obtained  throughout  are  applicable  also  for  sequential  launching  of 
narrow  band  pulses  Instead  of  postreception  frequency  splitting  of  a  wide 
band  pulse.  It  Is  known  (see  for  example  -  Beasely  and  Ward  (7))  that 
the  signals  resulting  from  this  procedure  are  statistically  independent. 
Since  the  bandwldths  of  the  signals  produced  In  this  way  are  n  times 
narrower  than  those  of  the  original  signal,  the  resulting  range  cells 
will  be  n  times  longer.  Hence  the  noise  power,  both  inside  and  outside 
the  target  range  cell,  will  be  n  times  larger  than  for  the  case  of  the 
Input  echo.  (It  is  shown  in  appendix  A,  however,  that  the  target  peak 
echo  location  is  unaffected  by  the  frequency  splitting  process  with  cer¬ 
tain  restrictions  on  the  spectra  of  the  transmitted  signal.  In  the 
appendix  we  use  Gaussian  shaped  signals.  However  this  restriction  that 
Is  Introduced  to  simplify  the  mathematics  Is  not  essential  as  In  success- 


^  Note  that  the  letter  n  when  subscripted  stands  for  noise  whereas  with¬ 
out  subscript  It  represents  the  number  of  frequency  splitting  channels. 
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ful  experiments  we  used  transducers  whose 

impulse 

response  envelope 

deviated  considerably 

from  Gaussian.) 

Thus 

our  two 

conditions 

to  be 

distinguished  for  the 

signals  resulting 

f  rr>E 

frequency 

splitting 

can  be 

written  as: 

r^  -  mfn*j|^ 

1  »  1 , . .  .n 

Ho=  ^i  -  “oi 

1  —  l,...n 

(3.6) 

where  rj  corresponds  to  the  received  signal  in  the  ith  window  at  a  certain 
delay,  m  is  the  target  amplitude  that  does  not  change  between  windows  and 
n*  and  n*  are  the  noise  components  corresponding  to  the  cases  (a  target 
in  the  range  cell)  and  (no  target  in  the  range  cell) .  Both  n*  and  n* 
are  Guassian  but  with  different  variances. 


The  variances  of  nj^^  and  n*j^  are  given  respectively  by 

*2  2 
®1  “  "«1 


After  squaring  each  of  the  channels,  the  SNR  of  any  one  channel  can 
be  written 

2  2- 

(Sj  «  ®  +  "'^1 

^  ICH  na^^  /3  (3.7) 

Summing  the  frequency  windows  and  averaging  we  obtain: 


^Lt  *^11  *^1^1  *  E  {  I  (nH-nJ^)^}  -  n(m2+aj^) 


1*1 


i*l 


=  n  /3n  0^2 

1-1 


(3.8) 

(3.9) 


and  the  output  signal  to  noise  ratio  becomes 
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2  ^  2 
+  nOj^ 


(3.10) 


Combining  eqts.  (5)  and  (10)  we  obtain  the  signal  to  noise  ratio  enhance¬ 
ment  (SNRE)  for  non-linear  averaging 


SNRE 


(S\  2  2 

''^‘'o  m  +  noj^ 

/n(m^  +  a?) 

jj  in  i 


/i-[(“  )^  i+  1] 

1  n _ 

(®_)^  +  1 
ni 


(3.11) 


Notice  that  this  enchancement  depends  only  on  the  clutter  noise 
variance  in  the  target  region,  and  that  it  varies  from  /n  at  m«a^  to 
l//n  at  m»oi.  Thus  signal  to  noise  ratio  enhancement  only  occurs  at 
small  input  SNR.  There  is,  however,  signal  to  noise  ratio  enhancement 
for  the  output  signal  compared  to  each  channel  of  the  split  signal  as  can 
be  seen  from  eqts.  (7)  and  (10)  which  give  for  this  case: 


SNRE  =>  /n 


(3.12) 


It  is  this  property  of  the  non-linear  averaging  algorithm  which  is  used  in 
frequency  agility  radar. 


Signal-to-Noise  Ratio  Enhancement  for  Linear  Averaging 

We  now  analyze  the  second  of  the  "classical"  frequency  compounding 
techniques  in  which  the  signals  produced  by  frequency  splitting  the  echo 
are  first  linearly  averaged  and  then  squared. 

It  is  obvious  that  splitting  the  signal  and  summing  the  results  will 
merely  restore  the  original  signal.  Thus  this  process  does  not  change 
the  relative  amplitudes  of  the  signal  with  respect  to  clutter.  Hence, 
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the  SNRE  of  linear  averaging  does  not  produce  any  Improvement  with  respect 
to  the  original  broad  band  echo.  The  algorithm  does,  however,  produce 
enhancement  with  respect  to  the  Individual  signals  produced  by  splitting 
the  echo  spectrum.  Thus  one  can  show  that  the  signal  to  noise  ratio 
Improvement  of  the  linearly  averaged  and  then  squared  signal  with  respect 
to  the  squared  signal  of  a  single  channel  Is 


n  y 


(3.13) 


The  fact  that  this  signal  to  noise  ratio  enhancement  approaches  n  at 
large  Input  SNR  explains  why  linear  averaging  Is  preferred  over  non-llnear 
averaging  In  frequency  agility  radar  when  the  Input  of  the  system  consists 
of  a  set  of  n  frequencies  from  which  a  wide  band  signal  can  be  constructed 
by  linear  summation. 


Slgnal-to-Nolse  Ratio  Enhancement  for  Minimization  of  Squared  Signals 

In  this  section  we  develop  the  statistics  of  the  minimized  signal  and 
show  that  the  Minimization  algorithm  provides  a  much  larger  Improvement  In 
signal  to  noise  ratio  than  the  two  techniques  analyzed  previously.  We 
start  as  before  by  splitting  the  received  echo  into  n  non-overlapping 
frequency  bands.  The  time  dependent  signal  at  each  of  these  frequencies 
Is  then  squared*  and  the  minimum  of  the  squared  signals  at  each  delay  Is 
chosen,  l.e. 


*In  order  to  allow  a  fair  comparison  between  the  techniques  analyzed  In 
this  paper,  we  analyze  here  the  effect  of  minimizing  the  square  of  the 
signals  produced  by  splitting  the  echo.  This  also  conforms  to  the  pro¬ 
cedure  used  for  the  simulated  and  experimental  data  which  we  compare 
with  our  analysis. 
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Z(t)  -  mlntr^(t) . 

In  order  to  Investigate  the  statistical  properties  of  the  random 
variable  z  -  mln[Tx...YQ]  where  Yi  represents  the  amplitude  of  the  squared 
signal  after  a  certain  time  delay  for  the  1th  frequency  window,  we  first 
construct  a  simple  experiment  In  which  we  choose  the  minimum  of  2  random 
variables 

Z-mln[Yi,Y2]  (3.14) 

It  Is  then  easy  to  show  (see  for  Instance  [2])  that  the  new  distribution 
function  of  z  Is  given  by 

PgC*)  “  ~  (3.15) 

where 

FY(y)  -  P(Y  <  y) 
and 

*'Yj^Y2(yi»y2)  “  F(Yl  <  71  ao*!  ^  72) 

If  we  create  a  random  variable  which  Is  the  minimum  of  n  random  variable 


Z  -  min  [Yl . Yn] 

Then  It  can  be  shown  by  Induction  that  the  probability  distribution  func¬ 
tion  F2(z)  will  be 


n  n-1  n  n-2  n-1  n 

Fz(z)  "  I  7y  (z)  -II  Fy.y/z.z)  +  I  I  I  Fy.Y  X  (z.a,z)- 

1-1  1  1-1  j-1+1  1  J  1-1  j-1+1  k-j+1  1  J  (3.16) 
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In  order  to  simplify  the  mathematics  we  will  assume  the  Y^'s  to  be  In¬ 
dependent  and  Identically  distributed  which  corresponds  to  no  overlap 
between  the  frequency  windows.  In  this  case  eq.  15  becomes 


F-(z)  -  1  -  (1  -  Fv(z))“ 


and  the  density  function  of  the  minimum  becomes 


(3.17) 


(z) 


nfY(z)(l  -  Fy(z)]“-1 


(3.18) 


If  we  split  the  received  signal  echo  Into  n  windows,  each  of  the  n 
time  signals  can  be  considered  to  be  a  Gaussian  process. 

After  squaring  the  signal  containing  clutter  only  the  density  func¬ 
tion  of  the  squared  amplitudes  for  each  window  becomes 


f§(y) 


-  y 

_ e  2^2 

/2iry'  0$ 


y  >  0 


(3.19) 


and  the  distribution  function  becomes 


FY(y) 


where 


2  erf  /X 
0%' 


y  >  0 


erf(y)  -  — L  /  e 

o 


2 

X 

2  dx 


(3.20) 


(3.21) 


The  minimized  clutter  noise  variance  can  now  be  evaluated  through 


a  2  »  J  z2  _  jjz  f£(z)dz]‘ 


(3.22) 


where  Is  the  density  function  of  the  minimized  squared  clutter 

signal  [eq.  18]  with  Fy(z)  ••  F^(z)  and  fY(z)  *  1^(2). 


When  a  target  exists  In  the  range  cell  the  statistics  of  the  returned 
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echo  changes  and  the  signal  can  no  longer  be  considered  as  zero  mean.  In 
addition  the  variance  of  the  noise  may  change  If  the  target  occupies  a 
significant  part  of  the  range  cell  thus  reducing  clutter  from  this  region. 
The  effect  on  the  target  echo  of  passing  It  through  different  frequency 
windows.  Is  analyzed  In  the  appendix  and  shown  schematically  In  Figure  2 
In  which  we  plot  the  target  echo  at  the  output  of  the  different  frequency 
windows.  One  can  see  that  the  maximum  of  the  target  signal  Is  unchanged 
In  height  and  occurs  at  the  same  time  delay  for  all  the  frequency  windows 
and  will  therefore  not  disappear  due  to  Minimization.  Hence  the  height 
of  the  target  signal  peak  Is  unchanged  by  minimizing  but  Its  width  Is 
reduced.  The  wldt^  Is  reduced  fflalnly  due  to  the  fact  that  away  from  the 
peak  region  the  signals  corresponding  to  the  different  windows  are  phase 
Incoherent  and  for  each  range  delay  away  from  the  peak  It  Is  highly 
probable  that  one  of  these  signals  will  be  very  close  to  zero.  The 
Minimization  which  chooses  the  smallest  signal  for  each  range  delay,  will 
therefore  narrow  target  echoes. 

The  received  target  signal  after  filtering  and  before  squaring  will 
have  the  following  density  function  for  a  target  in  the  range  cell 


1 

a*  e 


(x-m)' 


(3.23) 


where  Is  the  s.d.  of  the  clutter  In  the  range  cell.  The  density  func¬ 
tion  of  the  squared  signal  becomes  (for  a  target  In  the  range  cell) 


ft(y)  -  _.L_.  - 

2  /2Tiy  a*  [e  2a*2  g 


m)2 

2  0*2 


1 


y  >  0 

(3.24) 


the  distribution  function  becomes 


Target  Echo 
Amplitude 


B 


output 


Figure  3.2(A)  Target  echo  with  Gaussian  envelope  before  frequency  split¬ 
ting.  (B)  Target  echo  after  Hlnlmlaztlon,  showing  Invariance  of  peak 
position. 


the  mean  signal  of  the  target  can  be  evaluated  now  through 
00 

^out  "  /  (3.26) 

o 

where  as  In  eq.  (17) 

n-1 

f|(z)  -  nf5(z)  [1-f5(z)1  (3.27) 

likewise  the  output  noise  power  can  be  calculated  through 

— * 

“out^  "  /  (3.28) 

o 

where 

n-1 

f|(z)  -  nf|(z)(l-FC(z)J 

and  the  signal  to  noise  ratio  of  the  minimized  signal  is 


(3.29) 


Note  that  is  not  merely  a  function  of  the  target  signal  but  also  con¬ 
tains  a  contribution  due  to  the  squared  noise. 

Using  eqts.  (26)  and  (28)  the  signal  to  noise  ratio  enhancement  of 
the  Minimization  algorithm  has  been  calculated  numerically  and  is  shown  in 
Figure  3  as  a  function  of  the  input  SNR  for  the  case  of  oi  -  Oq  i.e.  with 
equal  noise  near  the  target  and  far  away  from  it.  Also  shown  in  Figure  3 
are  the  signal-to-nolse  ratio  enhancement  provided  by  the  two  conventional 
algorithms  as  described  earlier.  It  is  clear  from  the  figure  that  the 
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SNRE  for  linear  averaging  is  only  unity,  and  that  the  SNRE  for  the  non¬ 
linear  averaging  algorithm  is  even  worse. 

The  Minimization  algorithm  is  seen  to  be  highly  nonlinear.  With  high 
input,  SNR  the  output  SNR  is  significantly  higher  than  for  either  of  the 
conventional  techniques.  With  low  input  SNR  the  Improvement  is  not  as 
strong  and  its  performance  actually  falls  below  that  of  the  other  tech¬ 
niques  . 

Receiver  Operating  Characteristics 

The  signal  to  noise  ratio  enhancement  performance  is  undoubtly  an 
advantage  of  the  Minimization  algorithm  especially  if  we  intend  to  display 
the  processed  signal  without  further  processing.  In  this  case  the  ratio 
between  the  average  signal  to  the  noise  power  might  be  taken  as  a  crite¬ 
rion  of  the  visual  target  enhancement.  However,  if  we  are  interested  in 
target  detection  we  have  to  examine  the  ROC  of  the  proposed  algorithms. 
The  fact  that  the  signal  to  noise  ratio  obtained  by  the  Minimization 
algorithm  is  greater  than  those  obtained  by  the  other  algorithms  examined 
does  not  Insure  that  it  will  have  a  better  ROC.  This  is  illustrated  by 
Figure  4  which  shows  the  probability  density  function  of  the  output  of  a 
possibly  optimum  receiver,  compared  to  the  output  probability  density 
function  for  Minimization.  The  density  function  obtained  by  using  the 
Minimization  can  be  narrower,  providing  a  lower  variance  (lower  noise 
power).  However,  the  tail  of  the  probability  density  function  for  the 
minimized  signal  is  broader  than  that  of  the  optimum  receiver,  implying  a 
higher  probability  of  high  peaks  and  with  it  a  higher  probability  of  false 
alarm  for  Minimization  than  for  the  optimum  receiver.  We  assume  for 
illustrative  reasons  that  the  target  echo  size  was  not  significantly  re¬ 
duced  (Implying  a  high  signal  to  noise  ratio  in  the  input)  and  the  signal 
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Figure 

process 


3. A  Possible  probability  density  unctions  for  the  Minimization 
compared  to  that  of  an  optimum  receiver. 


52 


CO  noise  ratio  at  the  output  Is  determined  by  Che  amount  of  reduction  In 
the  noise  variance.  The  reduction  In  the  peak  target  echo  obviously 
effects  Che  detection  properties  adversely. 

To  calculate  Che  ROC  for  Minimization  we  use  the  model  of  eq.  6,  In¬ 
volving  uncorrelaCed  frequency  windows.  In  order  to  obtain  the  ROC  for 
the  correlated  case,  one  could  follow  the  approach  outlined  In  [42]. 
However,  as  stated  earlier  we  restrict  ourselves  here  to  the  case  where 
Che  different  frequency  windows  are  not  correlated.  It  Is  clear  chat 
squaring  Che  different  frequency  windows  will  not  change  Che  ROC.  Thus  It 
Is  sufficient  to  refer  to  the  absolute  value  of  the  signal  for  calcula¬ 
tions.  It  may  be  seen  Chat  for  a  given  threshold  z  the  probability  of 
detection  and  false  alarm  for  a  single  frequency  window,  becomes 

2  2 

-  ~  °)  -  (z  ~  °)  (3.30) 

•  2o  — z  2o 

P_(Di)  -  — L_  t/e  ^dz+  /e  ^  dz]«i  +  erf 

/2Tr  z  -«  ^  a* 

2 

_  z 

gs  20*^ 

P  (F,)  -  -Jl.  .  /  e  o  dz  -  1  -  2erf(£_  )  (3.31) 

where  Pr(Di)  Is  the  probability  of  detection  for  a  single  frequency  window 
and  Pr(^l)  ^3  the  probability  of  false  alarm  for  a  single  frequency  window 
for  a  given  threshold  z.  The  probability  of  detection  for  the  Minimiza¬ 
tion  process  equals  the  probability  that  the  smallest  of  n  target  echoes 
exceeds  the  threshold,  i.e.  the  probability  that  all  the  n  echoes  are 
equal  or  greater  than  z.  Hence  for  n  independent  windows,  the  probability 
of  detection  and  false  alarm  for  one  range  delay,  are  respectively, 

P^(D)  -  [Pr(Di)]“  -  [1  +  erf 

2  ol* 


(3.32) 
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slmlllarly,  the  probability  of  false  alarm  for  Minimization  is  equal  to 
the  probability  that  Independent  noise  signals  all  exceed  the  threshold. 
Thus 

n  n 

P-(F)  =  [Pr(Fi)]  =  [1  -  2erf  (3.33) 

Note  that  in  general  the  probability  of  false  alarm  Pp  of  the  whole  A-line 
depends  ou  the  exact  nature  of  the  fluctuations  of  the  clutter  scattering 
density  or  scattering  cross-section  and  should  be  defined  as 

00  00 

Pp  =  ^  Pr(F|ki)  P(ki)  =  I  Pr(F)P(ki) 

1=0  1=0 

where  kj  is  the  ratio  of  the  noise  std.  in  a  cell  without  a  target  to  the 
noise  std.  of  the  target  cell  and  PCk^)  is  the  probability  of  finding  such 
a  ratio  along  the  A-line.  P(kj^)  depends  on  the  properties  of  the  medium 
and  will  change  from  one  type  of  medium  to  another.  Thus  we  concentrate 
our  analysis  on  the  first  term  Py(F|kj^).  In  order  to  calculate  Pp  for  a 
specific  case  we  require  a  knowledge  of  P(kj^)  as  well  as  of  the  function 
Pj.(F|ki)  which  are  given  in  Figures  5(A)-(C)  below. 

The  ROC  for  Minimization  obtained  in  eqs.  32,  33  are  plotted  in 

Figure  5((A) , (B) , (C))  for  two,  four  and  six  windows  respectively.  The 
ratio  R-Oq/oi  of  noise  away  from  the  target  to  noise  near  the  target,  is 
used  as  a  parameter  in  these  curves.  For  the  sake  of  comparison  the  ROC 
calculated  for  linear  averaging  are  also  shown  (Averaging  is  known  to 
provide  the  best  ROC  for  the  case  Oq  =  oi).  One  can  see  that  the  ROC  for 
Minimization  Improves  with  respect  to  that  for  averaging  as  the  ratio 
between  Oq/oi  Increases.  For  example.  Figure  5(A)  shows  that  for  two  win¬ 
dows  averaging  is  preferable  to  Minimization  for  R<2  whereas  Minimization 
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Figure  3.5  Receiver  Operating  Characteristics  for  Minimization  and 
averaging  using  two,  four  and  six  frequency  windows  for  various  ratios 
of  Oq  to  01.  01  Is  the  s.d  at  the  target  range  delay  and  Oq  la  the  clutter 
s.d  elsewhere  ^curves  thlnned-out  In  the  top  right-hand  corner  to  retain 
clarity),  m/oj^  was  set  to  1. 

(A)  Two  frequency  windows. 

(B)  Four  frequency  windows. 

(C)  Six  frequency  windows. 
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Is  preferable  for  R>4.  Figure  4(B)  shows  Chat  Che  same  slCuacloa  holds 
for  four  windows.  Figure  5(C)  shows  that  for  n«6  Minimization  is  prefer¬ 
able  only  for  clutter  ratios  away  from  the  target  to  that  at  the  target 
of  R  >  8.  We  may  therefore  conclude  chat  if  Che  ratio  of  clutter  at  the 
target  to  clutter  elsewhere  is  small  enough.  Minimization  will  be  prefer¬ 
able  to  averaging  for  detection. 


EXPERIMENT 

The  theoretical  SNRE  plots  of  Figure  3  which  compare  Minimization 
with  linear  and  non-linear  averaging  are  found  to  agree  well  with  experi¬ 
mental  flaw-to-graln  echo  enhancement  curves  obtained  experimentally  in 
(4).  Flaw-to-grain  echo  ratio  is  defined  as  the  ratio  between  the  flaw 
(target)  peak  echo  amplitude  and  Che  largest  echo  amplitude  from  the 
clutter  present  either  in  the  squared  backscattered  echo  signal  (i.e. 
input  flaw-to-graln  echo  ratio)  or  the  processed  data  (i.e.  output  flaw- 
to-graln  echo  ratio).  Therefore,  this  is  an  analogous  but  not  quite 
identical  measurement  to  the  SNR  definitions  defined  for  the  three  algo¬ 
rithms  defined  above. 

The  experimental  data  shown  in  Figure  6  corresponds  to  a  signal  which 
is  obtained  by  summing  two  separate  signals:  1)  Grain  signals  obtained 
from  a  stainless  steel  sample  with  86  pm  average  grain  size,  and  ii)  A 
single  echo  from  a  flat  surface  reflector  simulating  the  flaw.  The  plots 
in  Figure  6  are  obtained  by  varying  the  amplitude  of  the  simulated  flaw 
signal  and  measuring  the  flaw-to-grain  echo  ratios  of  the  squared  received 
signal  and  the  processed  data  for  all  three  algorithms. 

Although  Che  particular  processing  parameters  selected  here  are  not 
necessarily  optimal  for  any  of  the  processing  algorithms,  they  result  in 
data  representative  of  the  general  behavior  of  each  algorithm  in  enhancing 
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the  flaw  echo.  A  more  detailed  descrlptloo  of  the  experimental  techniques 
used  to  obtain  Figure  6  Is  presented  In  (3,4). 

Bllgutay  In  his  experiments  used  at  Purdue  a  noise  correlation  system 
as  a  transmitter/receiver.  At  Drexel,  however,  we  used  a  pulser  as  a 
transmitter  and  digitized  the  received  echo  with  a  high  speed  digitizer. 
To  enhance  the  slgnal-to-thermal-nolse  ratio  we  averaged  the  received  echo 
by  summing  50  consecutive  returned  A  line  echoes.  Though  the  noise  cor¬ 
relation  system  output  signal  from  a  single  reflector  (see  (4))  Is  dif¬ 
ferent  from  the  received  signal  due  to  a  point  reflector  Illuminated  by 
a  pulse,  the  results,  however,  are  very  similar.  Experimental  results 
using  the  new  system  are  shown  In  Figure  7.  In  Figure  7(A)  we  combined 
the  received  echo  from  a  stainless  steel  with  average  grain  size  of  about 
76  ym  and  the  echo  reflected  from  a  flat  surface.  In  Figure  7(B)  the 
signal  of  7A  Is  squared  (to  enable  comparison  with  the  processed  signal 
which  Is  done  on  the  squared  windows)  and  In  Figure  7(C)  we  can  see  the 
signal  of  Figure  7(A)  processed  by  the  Minimization  Algorithm.  It  Is 
easy  to  see  that  there  Is  an  apparent  Improvement  In  a  slgnal-to-clutter- 
nolse  ratio.  We  did  not  attempt  here  to  optimize  the  slgnal-to-nolse 
ratio  by  changing  the  number  of  windows  or  the  amount  of  overlap  between 
them  as  was  done  for  example.  In  (4),  as  we  Investigate  here  a  theoreti¬ 
cal  model  and  not  a  particular  case.  Thus  the  signal  to  noise  ratio 
enhancement  In  our  case  Is  less.  The  advantage  of  using  the  new  system 
i.e.  using  a  pulser  and  a  high  speed  digitizer  Is  In  the  data  acquisition 
time.  The  noise  correlation  system  uses  a  water  delay  line  that  employs 
a  motor  to  control  the  delay.  As  a  result  acquisition  of  a  single  A  line 
composed  of  512  sample  points  takes  about  20-30  seconds.  If  the  tech¬ 
niques  Introduced  In  this  research  work  are  to  be  Industrially  Implemented 
a  high  speed  acquisition  system  Is  clearly  a  must. 
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Figure  3.7(A)  A  signal  resulting  from  the  sum  of  a  signal  from  a  stainless 
steel  (with  an  average  grain  size  of  75  um)  and  a  flat  surface  reflector 
target . 

(B)  The  signal  of  Figure  3.6A  squared. 

(C)  Processed  output  for  the  Minimization  Algorithm  for  10  non-overlapping 
frequency  windows. 
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It  is  clear  from  Figure  6  that  for  both  averaging  algorithms  the  out¬ 
put  flaw-to-grain  echo  ratio  is  linearly  dependent  on  the  input  flaw-to- 
graln  echo  ratio,  which  indicates  that  the  performance  of  these  algorithms 
is  not  strongly  dependent  on  the  variation  in  flaw  amplitude.  These  re¬ 
sults  also  indicate  that  the  linear  averaging  algorithm  performs  slightly 
better  than  the  non-linear  averaging  algorithm,  which  agrees  with  the 
theoretical  results  shown  in  Figure  3.  However,  neither  averaging  algo¬ 
rithm  shows  any  noticeable  enhancement  in  flaw-to-graln  echo  ratio  experi¬ 
mentally,  in  agreement  with  the  theoretical  derivations. 

The  Minimization  algorithm  is  seen  to  show  more  sensitivity  to  the 
flaw  amplitude  (as  is  evident  from  the  non-linear  nature  of  the  input- 
output  flaw-to-graln  echo  ratio)  than  the  averaging  algorithms.  These 
results  are  also  in  close  agreement  with  the  theoretical  plot  shown  in 
Figure  3. 

DISCUSSION 

In  this  chapter  we  analyzed  and  compared  the  performance  of  three 
frequency  compounding  algorithms  in  Improving  the  ultrasound  visibility 
and  detectability  of  targets  hidden  by  additive  clutter  or  speckle,  using 
only  the  Information  presented  in  a  single  A-scan. 

Many  criteria  need  to  be  considered  when  comparing  algorithms  for 
improving  target  detectability.  These  include  signal-to-nolse  ratio  or 
contrast,  probability  of  detection  and  false  alarm  (ROC),  bias  of  signal 
estimates,  dynamic  range,  longitudinal  and  lateral  resolution,  and  sim¬ 
plicity  of  implementation.  Each  of  these  criteria  can  be  improved,  but 
often  only  at  the  expense  of  increased  dynamic  range  requirements,  and 
the  Introduction  of  bias  in  target  amplitude  estimation. 
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In  this  chapter  we  analyzed  the  first  two  of  these  criteria  for  the 
three  algorithms  under  consideration,  namely  slgnal-to-nolse  ratio  and 
receiver  operating  characteristic.  It  was  shown  theoretically  and  con¬ 
firmed  by  experiment,  that  the  slgnal-to-nolse  ratio  or  contrast  enhance¬ 
ment  produced  by  the  Minimization  algorithm  Is  far  greater  than  that  of 
either  linear  or  nonlinear  averaging,  provided  that  the  Input  slgnal-to- 
nolse  ratio  Is  larger  than  unity.  This  holds  true  whatever  the  ratio 
between  the  clutter  density  In  the  target  region  and  elsewhere.  This 
advantage  of  Minimization  with  respect  to  linear  or  nonlinear  averaging  la 
however,  to  some  extent  Illusory,  since  It  may  be  accompanied  by  Increased 
probabilities  of  "target  drop-out".  In  other  words  even  though  the  con¬ 
trast  and  visibility  of  the  target  Is  Improved  on  the  average  by  Minimiza¬ 
tion,  there  la  also  an  Increased  probability  that  the  target  echo  will  be 
supreased  by  the  subtraction  of  an  awkwardly  located  noise  spike. 

These  consideration  are  quantitated  by  the  calculated  receiver  oper¬ 
ating  characteristics  Illustrated  In  Figure  5  which  show  that  for  a  given 
false  alarm  rate,  the  probability  of  detection  for  Minimization  Is  only 
Improved  over  that  for  linear  averaging,  when  the  clutter  density  near 
the  target  Is  at  least  four  times  smaller  than  that  elsewhere.  This 
situation  Is  In  fact  not  uncommon,  since  In  attenuatlve  media,  the  clutter 
echo  from  regions  close  to  the  transducer  Is  larger  than  that  from  deeper 
lying  target  regions  and  so  might  be  the  case  In  fluctuating  random  media. 
This  situation  Is  not  rectified  by  the  use  of  time-varying  gain  except  In 
the  relatively  rare  cases  where  the  attenuation  and  clutter  density  are 
known  at  every  point  of  the  medium. 

It  lo  111  luct  Ignorance  of  the  local  properties  of  the  medium  which 
limits  the  effectiveness  of  our  procedures.  In  the  forthcoming  chapter  we 
derive  the  optimum  receiver  for  the  target  model  analyzed  In  this  paper, 
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RESOLUTION:  Since  Minimization  Involves  filtering  the  target  echo  with  a 
narrow  bandwidth  filter,  it  reduces  the  longitudinal  resolution  in  the 
sense  that  target  echoes  may  enhance  or  depress  each  other  over  a  distance 
c/2AF  where  c  is  the  velocity  of  sound  and  AF  the  filter  bandwidth.  Thus 
a  strong  target  might  partially  suppress  an  adjacent  smaller  target. 

However,  because  of  the  effect  discussed  in  connection  with  Figure  2, 
Minimization  narrows  the  processed  target  echo  to  about  one  wavelength  of 
of  the  highest  center  frequency  for  splitting  the  spectrum,  i.e.  produces 
sharpening  in  the  processed  image. 


CONCLUSIONS 

In  this  chapter  we  analyzed  the  so-called  Minimization  algorithm  and 
compared  its  performance  in  Improving  the  visibility  of  targets  hidden 
by  additive  time-invariant  clutter  or  speckle  noise,  with  that  of  two 
"classical"  algorithms  used  in  frequency  agility  radar,  linear  and  non¬ 
linear  averaging.  All  three  algorithms  "frequency  compound"  signals 
produced  by  filtering  a  single  A-scan  ultrasound  echo  with  a  number  of 
non-overlapping  spectral  "windows".  The  operation  of  these  frequency 
compounding  algorithms  is  therefore  relatively  simple  since  they  operate 
on  only  one  A-mode  echo  at  a  time. 

Minimization  was  shown  to  be  capable  of  producing  significant  improve¬ 
ment  in  SNR  and  is  predicted  to  produce  sharpening  in  2-D  images.  Neither 
of  these  effects  are  produced  by  the  other  two  algorithms  studied. 

Linear  averaging  la  known  to  be  the  optimum  detection  algorithm  for 
additive  noise  and  accordingly  was  found  to  give  the  beat  ROC  of  the  three 
algorithms  tested.  However,  Minimization  is  shown  to  produce  better  ROC 
when  the  clutter  near  the  target  has  a  standard  deviation  which  is  at 
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least  four  times  smaller  than  that  elsewhere.  Such  cases  are  often  en¬ 
countered  In  attenuating  media  due  to  shadowing  by  large  targets,  and  for 
media  which  exhibit  fluctuations  In  either  scattering  concentration  or 
scattering  cross-section  along  the  A-llre.  The  analyses  In  this  chapter 
have  been  restricted  to  non-overlapping  frequency  windows.  The  conclusions 
reached  should,  however,  also  apply  to  filtering  with  overlapping  windows, 
as  Is  shown  experimentally  in  (18)  and  will  be  further  discussed  In 
chapter  5. 


CHAPTER  4 


OPTIMAL  PROCESSING 


INTRODUCTION 

Chapter  3  was  dedicated  to  the  analysis  of  the  minimization  algorithm. 
We  found  out  that  this  algorithm  provides  a  significant  Improvement  In 
signal  to  noise  ratio  with  a  marginal  loss  In  dynamic  range.  However,  the 
Receiver  Operating  Characteristic  (ROC)  was  found  to  be  less  Impressive. 
Only  for  a  high  ratio  of  variances  between  a  cell  containing  clutter  and 
the  target  range  cell  was  the  ROC  found  to  be  superior  to  averaging. 

In  this  chapter  we  Introduce  a  new  technique  for  Improving  the  ROC 
of  the  received  A-llne.  This  technique  which  we  call  "optimal  processing" 
Is  also  based  on  split  spectrum  processing.  In  the  first  part  of  this 
chapter  we  Introduce  the  theoretical  basis  and  motivation  for  split 
spectrum  processing.  We  find  that  In  order  to  construct  the  optimum 
receiver  (a  receiver  that  will  minimize  the  probability  of  error  for  a 
given  threshold)  we  should  split  the  spectrum  Into  Independent  frequency 
windows  whose  number  Is  directly  proportional  to  the  required  range 
resolution.  From  the  set  of  random  variables  that  correspond  to  each 
range  delay  we  construct  the  receiver  algorithm.  For  trivial  cases  i.e. 
for  cases  where  the  noise  and  signal  spectra  are  white  In  the  frequency 
range  of  Interest,  the  algorithm  becomes  the  estimated  mean  divided  by 
the  estimated  standard  deviation  of  the  spectral  decomposition  components. 

We  proceed  to  evaluate  the  ROC  of  the  new  algorithm  for  media  with 
fluctuating  clutter  noise  and  compare  It  with  the  ROC  obtained  by  thresh- 


oldlag  the  received  signal  without  special  processing.  Then  we  introduce 
a  simplified  version  of  the  mathematical  model  for  which  we  replace  the 
matched  filters  simrT(fQ-f )/ir(f ji~f )  obtained  in  the  optimum  receiver 
derivation,  with  Gaussian  shaped  filters  and  show  that  its  performance 
is  very  close  to  the  performance  of  the  mathematical  model.  We  then 
implement  this  simplified  algorithm  and  show  experimentally  its  effective¬ 
ness  in  target  enhancement.  Finally  we  evaluate  experimentally  the 
resolution  performance  of  the  new  algorithm  in  comparison  to  Minimization. 

THEORY 

As  is  known,  the  received  clutter  noise  power  varies  as  a  function 
of  range  delay  due  to  changes  In  scattering  density,  scatte^^ng  cross- 
section,  attenuation,  beam  Intensity  etc.  However  if  we  take  a  small 
enough  time  Interval  [Ti,Tf]  of  the  received  echo  e(t)  we  can  assume  in 
most  practical  cases  that  the  statistical  properties  of  the  echo  in  that 
interval  are  constant.  We  will  treat  the  problem  of  clutter  reduction 
in  this  Interval  [Tj:,Tf]  as  a  problem  in  statistical  decision  theory. 
The  general  Hypothesis  set  can  be  written  as  follows, 

H^:  r(t)  -  ux(t-Tt)  +  a(t)  T^  <  t  <  Tf  4.1 
Hq:  r(t)  -  n(t)  Tf  -  Ti  -  T 

where  n(t)  is  a  Gaussian  random  process  with  variance  a^,  a,  p  and  Tt  are 
unknown  and  x(t)  is  a  known  signal.  Though  the  returned  echo  [-«,  <*>]  is 
in  general  not  stationary,  we  assume  that  in  the  interval  of  analysis 
[Ti,Tf]  the  noise  process  n(t)  would  be  a  sample  function  of  a  stationary 
process  (a  requirement  that  la  met  if  the  Interval  [Ti,Tf]  is  sufficiently 
small).  The  problem  is  to  find  a  decision  rule  as  to  whether  p^O  or  p^O 
where  a,  p  and  tj-  are  unknown. 
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Let  us  consider  first  the  problem  for  it  known.  In  this  case  as  we 
shortly  show,  the  decision  rule,  can  be  constructed  from  the  set  of  random 
variables  obtained  by  an  orthonormal  expansion  (Karhunen-Loeve  expansion, 
see  for  example  (42))  of  the  random  process  r(t).  In  this  technique  we 
transform  the  problem  from  a  random  process  problem  into  a  random  vari¬ 
ables  problem  (see  for  example  Scharf  (43)).  These  random  variables  are 
the  r^’s  in  the  expression 

n 

r(t)  -  I  ri  4ii(t)  4.2 

i=0 

where  the  functions  i>i(t)  are  orthonormal.  To  explain  the  procedure  of 
obtaining  the  coefficients  of  the  expansion  (the  random  variables)  we 
assume  for  simplicity  that  the  spectrum  of  n(t)  is  rectangular  as  shown 
in  Figure  1,  and  that  the  spectrum  of  x(t)  is  no  larger  than  B  (There  is 
an  apparent  contradiction  between  the  time  limited  interval  of  the  ana¬ 
lyzed  signal  r(t)  and  the  limited  bandwidth  of  the  rectangular  noise 
process  assumed  in  the  example.  However  this  assumption  is  often  made 
to  simplify  an  analysis  and  in  most  cases  is  found  to  be  justified). 

To  find  the  '|)i(t)  one  solves  the  Integral  equation, 

T 

^i'}'i(t)  -  Rjjjj(t-T)  i)i^(T)dT  4.3 

2 

where  Is  both  the  eigenvalue  and  noise  power  associated  with  the  1-th 
orthonormal  function,  and  Rnti(T)  is  te  autocorrelation  function  of  the 
noise  process.  One  can  show  that 

t.(t)  =  i  e^  ^  for  f.  »  if  ;f»i  -Z<t<X  4.4 

^  /f  1  o  o  .j.  2  2 
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and  that 

Xi  s  Sn(fi)  4.5 


For  our  choice  of  spectrum  (shown  in  Figure  1),  all  the  Xj^'s  and 
thus  the  variances  of  all  the  orthonormal  coordinates  are  identical.  This 


makes  them  identically  distributed  (if  we  have  only  clutter  noise  in  the 
input).  Naturally  we  are  not  Interested  in  the  orthonormal  functions 
outside  the  frequency  band  B  of  the  clutter  noise  since  they  carry  no 
information.  (We  assume  that  the  spectrum  of  x(t)  is  no  broader  than 
that  of  n(t)  so  that  there  is  no  energy  outside  of  the  frequency  band  of 
width  B) .  Note  that  for  the  general  case  where  Sji(f)  is  not  white  in 
the  frequency  range  of  Interest  the  r^’s  will  not  be  identically  distri¬ 
buted.  The  assumption  of  a  white  spectrum  in  the  example  is  therefore 
merely  used  to  provide  an  insight  to  the  general  solution  of  the  problem. 


In  order  to  obtain  the  set  ri 's  of  the  random  variables  associated 
with  the  K.L  expansion  we  have  to  perform  the  process  shown  in  Figure  2. 
Since  the  number  of  filters  i^j^(T-t)  and  the  number  of  inde¬ 
pendent  random  variables  directly  proportional  to  the  observed 
time  T  and  equals, 


n  »  B  T 


4.6 


The  chosen  time  window  [Tj^.Tf]  determines  the  resolution  because  we  must 
assume  that  only  a  single  target  exists  in  this  window,  since  otherwise 
the  assumption  that  the  target  signal  x(t)  is  known  would  no  longer  be 
valid. 

Now  for  the  assumption  that  the  target  location  in  r(t)  is  known  we 
construct  the  likelihood  ratio  test  (LRT)  (The  LRT  gives  a  decision 
criterion  subject  to  a  certain  risk  function  to  determine  which  of  two 


Figure  4.2  GeneraCloa  of  expansion  coefficients, 

(A)  Correlation  operation. 

(B)  Filter  operation. 
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possible  hypotheses  is  true.  It  Is  composed  of  the  ratio  of  the  prob¬ 
ability  density  functions  of  the  observables  of  the  two  hypotheses.)  from 
the  set  of  random  variables  r^'s-  The  random  process  problem  set  in 
eq.  (1)  can  now  be  restated  as  a  random  variables  problem 

Hi:  ri  =  p  Xi  +  ni 

4.7 

Hq:  ri  =  ni 

where  the  x^'s  are  the  spectral  decomposition  components  r^  of  the  signal 
x(t)  fed  into  the  system  shown  in  Figure  2.  Instead  of  solving  directly 
the  problem  for  which  the  noise  spectrum  is  white  in  the  frequency  range 
of  Interest  (as  in  our  example)  we  quote  first  the  general  solution  to 
the  problem  of  colored  Gauslan  noise.  Later  we  return  to  our  example  as 
a  special  case  of  the  general  solution. 

If  the  spectrum  of  the  received  clutter  signal  is  not  white  in  the 
frequency  range  of  interest  or  the  signal  was  not  whitened  before  decom¬ 
position  (the  signal  can  be  also  whitened  during  decomposition  by  assign^ 
ing  appropriate  constant  multipliers  to  the  filters  (fri(t)),  the  likeli¬ 
hood  ratio  test  (LRT)  becomes  (see  for  example  Scharf  and  Lyte  (43), 
Helstrom  (15)), 


L(r) 


n 

I  xi^^i/^i 
1-1 


1 


n-1 


(I  ■  t 

1-1 


i-l  H 


4.8 


where  is  the  eigenvalue  associated  with  the  1-th  orthonormal  coordinate 
and  Xi  is  as  above.  The  statistics  of  L(0  is  the  student  t  distribution 
and  we  will  discuss  some  of  its  properties  later  in  the  chapter.  For  the 


r - - 

.  m 
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processing  described  here 

we  whiten  the  clutter  spectrum  making  It  similar 

’1 

m 

to  the  example  and  assume  that  the  target  energy  spectrum  Is  similar  to 

the  clutter  power  spectrum,  which  means  that  the  target  Is  assumed  to  be 

a  strong  point  reflector.  We  can  then  scale  x(t)  such  that  l  “  1  so 

1-1  H 

that 

n 

U  I  -  y 

1-1 

4.9 

i 

but  Xj-x  and  for  a  white  spectrum  within  the  frequency  band  B  so 

2 

or  -1 

X 

• 

\  n 

4.10 

1 

and  the  LRT  becomes 

1  r 

-  i.  *^1 

n  1.1 

L(r)  - 

<Jr 

4.11 

This  Is  nothing  else 

but  the  estimated  mean  of  the  spectral  decom- 

position  components  ri's 

divided  by  their  estimated  standard  deviation 

m 

and  can  be  written 

m 

L(r)  -  _L 
®r 

4.12 

m 

To  estimate  which 

Is  the  noise  power  after  whitening  we  have  to 

calculate 

i 

•^r  ”  t  ^  (  I  I  ^ 

n-1  1*1  n  1.1 

■ll')}''"  4.13 

.  « 
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and  the  LRT  becomes 
n 

^  l^i 

n  i«i 

L(r)  »  - 


4.14 


n-1  i=i 


I  r,r]] 

u 


2^,l/2 


As  the  number  n  of  the  orthonormal  functions  Increase  the  performance  of 
the  algorithm  In  estimating  a  Improves.  Note  that. 


n  n  Tf  Tf 

-  I  -  I  /  r(t)  (ti.(t)dt  =  i  /  r(t)  I  (|..(t)dt 
n  1-1  ^ 


4.15 


1=1  T, 


Recall  that  for  eq.  4.14  to  hold  ^  =  A 


n  n 

X(t)  “  I  Xi'l'iCt)  -  /A  ^  (tlj.(t) 
1*1  °  1=1 


so 


I  -  /-  x(t) 
1=1  ^ 


and  eq.  4.12  bec'!res 


n  _  Tf 

I  1^1  =  /  r(t)x(t)dt 

1=1  *  T, 


1=1 


4.16 


4.17 


4.18 


which  Is  seen  to  be  an  estimate  of  the  energy  of  x(t)  In  r(t).  So  eq.(14) 
Is  just  the  value  of  the  estimated  signal  energy  at  the  analyzed  range 
delay  divided  by  the  estimated  standard  deviation  of  the  clutter  in  [Tf, 
Tf].  Naturally  as  the  number  of  elements  In  the  expansion  increases  so 
does  the  accuracy  of  the  variance  estimate  and  the  performance  of  the  LRT. 
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Now  that  we  have  derived  an  optimum  LRT  to  decide  whether  an  individu¬ 
al  range  cell  does  or  does  not  contain  a  target,  we  return  to  the  original 
problem  in  which  we  treat  the  A-line  as  a  collection  of  m  range  cells. 
Each  of  these  corresponds  to  a  certain  range  delay  of  length  T  and  pro¬ 
duces  n  spectral  decomposition  components.  Since  we  assumed  n(t)  to  be 
Gaussian  these  n  r^  terms  for  each  cell  are  samples  from  a  Gaussian 
distribution.  It  is  important  to  note  that  the  cell  statistics  can  change 
from  cell  to  cell  due  to  either  a  change  in  scatterer  density  which  leads 
to  a  change  in  the  variance  of  r(t)  or  due  to  the  existence  of  a  target 
in  the  analyzed  time  window  which  leads  to  a  change  in  its  mean.  (In 
chapter  2  we  showed  that  such  a  change  (gradient)  should  lead  to  a  co¬ 
herent  reflection.  However  if  the  change  is  slow  compared  to  a  wave¬ 
length  (as  is  the  assumption  here)  the  coherent  component  is  negligible). 
Recall  that  we  do  not  have  any  Information  on  the  location  of  the  target. 
Thus  we  have  to  find  a  LRT  for  every  range.  There  are  three  possibilities, 

1.  The  analyzed  range  delay  is  such  that  T=Tt  where  Is  the  target 

range  delay.  This  case  is  referred  to  as  H^  and  was  analyzed  earlier. 

2.  There  is  no  target  in  the  time  window  of  length  T.  This  case  was 
also  analyzed  and  is  referred  to  as  Hq. 

3.  There  is  a  target  in  the  time  window  T  but  For  now  we  assume 

that  the  resolution  cell  is  no  smaller  than  the  analysed  time  Interval 
T.  In  this  case  multiple  targets  in  the  time  interval  T  should  be 
Interpreted  as  a  single  target.  We  will  return  to  the  resolution 
performance  of  this  algorithm  later  in  this  chapter. 

If  we  presented  the  likelihood  ratio  as  a  function  of  range  delay  we 

would  obtain  a  large  signal  wherever  a  target  exists  in  the  analyzed  range 

cell.  In  the  absence  of  a  target  the  estimated  mean  is  expected  to  be 
small  and  the  processed  output  would  consequently  be  small.  We  would 
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also  be  able  to  set  a  threshold  and  reject  signals  below  a  prespecified 
signal  to  noise  ratio. 

To  maximize  the  resolution  of  the  algorithm  (eq.  14),  set  the  fre¬ 
quency  windows  should  overlap  so  that  they  can  be  as  wide  as  possible. 
The  n  random  variables  corresponding  to  a  certain  range  delay  would  then 
become  correlated.  If  the  filters  are  correlated  then  the  correlation 
between  any  two  channels  containing  only  clutter  can  be  shown  to  be 
(assuming  that  the  signal  Is  whitened  before  processing  so  the  clutter 
power  In  all  the  channels  Is  Identical) 


r  Fi(f)Fj(f)dt 

■A  ^ 

r  lFi(f)l2df 


4.19 


For  Gaussian  windows  (see  (4)) 
f,-f,  2 

-  In  16(_ii) 

Pj  j  “  e  4.20 

where  b  Is  the  bandwidth  of  the  windows  and  f^,  f j  are  the  center  frequen¬ 
cies  of  windows  1  and  j.  The  LRT  for  this  case  would  be  (19) 

L(r)  =  xT  K"1  £  4.21 

where  K  Is  the  covariance  matrix  which  Is  composed  of  the  above  j's  and 
X  Ls  the  mean  vector  of  the  set.  This  can  be  obtained  by  feeding  the 
system  In  Figure  1  the  signal  x(t),  with  the  coefficients  ri,...rn  In  the 
figure  Interpreted  as  Xj^,...Xjj.  The  vector  f  Is  composed  of  the  set 
ri,...rn  obtained  by  feeding,  the  system  In  Figure  1  the  signal  r(t).  To 
construct  K  we  need  also  to  know  the  variance  of  the  channels 
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*1 

1  ^12  •• 

*2 

^12  1 

• 

*3 

K  -  0^ 

•  • 

r  - 

• 

_*IL_ 

_  .  .  1  _ 

_  ^^n  _ 

a  can  be  estimated  from  the  r^'s  following  the  same  approach 
as  before  (eq.  13),  If  the  noise  spectrum  Is  white.  For  colored  noise 
we  should  whiten  the  signal  before  processing  and  then  apply  eq.  14. 

It  Is  shown  later  that  the  range  resolution  of  the  optimum  receiver 
algorithm  Is  not  significantly  degraded  In  comparison  to  the  unprocessed 
signal,  even  when  the  frequency  windows  do  not  overlap.  We  will  therefore 
not  consider  the  case  of  overlapping  frequency  windows  for  target  enhance¬ 
ment  further. 

Detection  Performance 

The  next  step  Is  to  evaluate  the  performance  of  the  optimal  processor 
which  Is  optimum  In  the  sense  that  It  provides  the  best  probability  of 
detection  for  a  given  probability  of  false  alarm  for  the  analyzed  range 
cell  as  defined  In  eq.  (14). 

We  recall  that  In  our  procedure  the  original  signal  Is  segmented 
Into  m  cells  and  then  spectrally  decomposed  to  obtain  n  random  variables 
from  which  the  LRT  Is  calculated  as  depicted  In  Figure  2.  Now  let  us 
assume  a  target  somewhere  along  the  A-llne.  We  define  the  Probability 
of  Detection,  P(j  as  the  probability  that  the  processed  data  of  the  cell 
containing  the  target  will  exceed  a  certain  threshold.  The  probability 
cf  tc Ise  alarm  Pf  Is  defined  as, 
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where  P(f|ki)  is  the  probability  that  the  processed  data  of  the  i-th 
cell  having  a  variance  times  that  of  the  target  containing  cell  will 
exceed  the  threshold  set  for  the  k-th  cell  (the  target  cell),  p(k^)  is 
just  the  probability  of  finding  such  a  ratio  along  the  A-line.  Naturally 
for  each  scenario  of  fluctuations  appearing  along  the  A  line  a  different 
Pf  will  exist.  In  the  next  section  we  will  describe  the  detection  per¬ 
formance  of  the  optimal  processor  for  different  kj^-s. 

It  is  clear  that  the  performance  of  an  algorithm  depends  strongly 
on  the  distribution  of  p(ki^).  Naturally  we  cannot  construct  an  algorithm 
that  will  be  optimal  for  all  possible  fluctuations  in  the  clutter  variance 
along  the  A-llne  (i.e  p(ki)).  A  rough  measure  of  comparison  between 
algorithms  might  be  to  check  P(f]kj^)  for  the  worst  possible  performing 
cell  (the  cell  with  the  highest  variance)  and  probability  of  detection 
for  the  cell  with  the  lowest  clutter  power.  In  order  for  this  test  to 
be  meaningful  the  algorithm  should  be  monotonlc  in  the  sense  that  for  2 
cells  with  >  02  the  probability  of  false  alarm  of  the  processed  data 
of  cell  1  will  be  equal  or  greater  than  the  probability  of  false  alarm 
of  the  processed  data  In  cell  2. 

It  is  Important  to  note  that  the  best  detection  procedure  will  not 
necessarily  provide  the  best  Receiver  Operating  Characteristics  (ROC)  for 
this  problem.  (Recall  that  ROC  is  a  plot  of  the  probability  of  detection 
as  a  function  of  probability  of  false  alarm  of  a  receiver  for  varying 
thresholds  where  usually  the  signal  to  noise  ratio  Is  a  parameter).  We 
will  give  a  simple  example  for  this  argument.  If  all  the  m  cells  had  the 
same  variance,  the  optimum  algorithm  is  the  sum  of  the  elements  In  each 
of  the  cells.  Since  the  variance  in  all  the  cells  is  the  same,  the 
variance  in  the  LRT  becomes  a  proportionality  constant  and  the  LRT  reduces 
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to  the  sum  of  the  elements  in  each  of  the  cells .  The  variance  now  no 
longer  Influences  the  ROC,  it  merely  changes  the  threshold  for  which  the 
probability  of  detection  and  the  probability  of  false  alarm  is  computed. 
This  example  shows  that  the  best  detection  scheme  for  each  of  the  cells 
(Which  is  composed  of  the  estimated  mean  divided  by  the  estimated  standard 
deviation)  is,  in  general,  not  the  beat  algorithm  for  improving  the  ROC  of 
the  whole  A  line.  The  conclusion  is  that  our  optimum  receiver  has  to  be 
checked  against  other  algorithms  for  different  ratios  of  variance  in  the 
target  range  cell  and  the  variance  elsewhere.  Depending  on  the  actual 
conditions,  an  optimum  algorithm  can  then  be  chosen. 

I  *  we  knew  the  ratio  k  between  the  variances  in  the  target  range  cell 
and  the  cell  for  which  P(flkj^)  is  calculated,  then  an  optimum  receiver 
could  be  constructed.  To  derive  the  LRT  for  this  case  we  use  eq.  327  in 
chapter  2  vol.  1  of  "Detection,  Estimation,  and  Modulation  Theory"  (42). 
This  equation  is  the  general  binary  solution  for  the  LRT  for  Gaussian 
random  variables  with  different  vector  mean  and  covariance  matrices.  In 
our  case  the  random  variables  are  Independent  and  the  LRT  reduces  after 
some  algebra  to, 

n  n 

1  I  rjd-k)  +  li(  I  r^)2  -  L(r)  4.24 
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where  k  is  the  ratio  of  the  variance  of  the  cell  containing  clutter  only 
to  that  of  the  target  cell. 

In  this  case  the  probability  of  false  alarm  will  be  determined  by  the 
cell  containing  only  clutter  and  the  probability  of  detection  will  be 
determined  by  the  target  size  and  clutter  noise  in  the  target  range  cell. 
p(f|ki)  and  pj  are  calculated  through. 
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where  and  Lq  are  L(f)  for  and  Hq  respectively 

Eq.  24  shows  that  the  LRT  can  range  from  the  sum  of  the  r^  terms  in 
each  cell  for  k=*l  to  a  close  approximation  of  the  variance  of  r^  for 
k»l.  Naturally  choosing  a  given  k  for  constructing  the  receiver  would 
worsen  the  performance  for  any  ratio  which  is  greatly  different  from  the 
chosen  one. 

In  general,  the  probability  density  p(k£)  of  the  clutter  in  the 
various  cells  of  the  A  line  will  vary  from  one  A  line  to  the  next,  and 
will  usually  be  unknown.  Thus  there  is  usually  not  enough  information  to 
construct  an  optimum  receiver  for  a  whole  A  line,  and  even  if  there  were, 
the  form  of  the  algorithm  would  usually  have  to  be  changed  for  each  A 
line  which  would  be  an  extremely  impractical  procedure.  We  will  therefore 
in  most  cases  be  constrained  to  restrict  ourselves  to  the  LRT  of  eq.  12 
which  provides  the  maximum  probability  of  detection  for  a  given  prob¬ 
ability  of  false  alarms  for  a  target  in  one  single  range  cell. 

Before  analyzing  the  performance  of  the  LRT  we  point  out  an  important 
property  of  the  function  L(f)  which  is  identical  with  the  well  known 
Student  t  distribution.  For  this  distribution 
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where  fL  (L^)  is  the  density  function  of  L(£)  if  a  target  exists  in  the 
range  cell  and  fr  (L-)  is  the  density  function  of  L(?)  if  no  target  exists 
in  the  range  cell  (w“0). 

This  distribution  can  be  seen  to  have  the  property  that  the  prob¬ 
ability  density  function  (Lq)  is  only  a  function  of  n,  the  number  of 

o 

windows.  Thus  for  y-0,  i.e.  in  the  absence  of  a  target  the  variance  of 
L(0  is  Independent  of  the  variance  of  the  spectral  decomposition  com¬ 
ponents  rjL  •  One  of  the  consequences  of  this  property  is  that  even  though 
an  A  line  may  exhibit  strong  noise  fluctuations  in  different  range  cells, 
its  LRT  !(!?■)  will  have  constant  variance  in  every  range  cell.  Another 
consequence  is  that  using  eq.  12,  the  false  alarm  rate  for  cells  not 
containing  targets  for  a  given  threshold  is  Independent  of  the  noise  in 
these  cells.  This  is  clearly  a  very  useful  property. 

The  number  of  variables  in  the  LRT  equation  (eq.  11)  is  too  large 
to  fully  demonstrate  its  detection  properties  on  one  graph.  We  chose 
therefore  to  set  a  constant  false  alarm  probability  and  plot  the  prob¬ 
ability  of  detection  as  a  function  the  input  of  signal  to  noise  voltage 
ratio  in  the  target  range  cell.  False  alarm  probability  Pf  and  prob¬ 
ability  of  detection  are  determined  from  aqs.  25  and  26,  using  tables 
of  the  student  t  distribution.  The  parameter  was  chosen  to  be  the  number 
of  elements  in  the  expansion  (or  the  number  of  Independent  windows).  In 
Figure  3  and  4  we  present  the  theoretical  results  for  several  n  the  number 
of  windows.  It  is  clear  that  Increasing  the  number  of  windows  Improves 
the  ROC.  Of  course  as  n  Increases  the  width  of  each  spectral  window 
decreases  which  increases  the  time  window  [Ti^,Tf]  and  then  worsens  the 
resolution.  Figure  3  correspond  to  a  constant  probability  of  false  alarm 
(CFAR)  of  0.01  and  Figure  4  correspond  to  a  CFAR  of  0.05.  One  can  im¬ 
mediately  observe  that  as  the  number  of  elements  n  in  the  expansion 
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Figure  4.3  Detection  probability  versus  slgnal-to-nolse  voltage  ratio 
for  optimal  processing  for  various  expansion  sizes  n.  The  probability 
of  false  alarm  la  set  to  0.01.  (taken  from  36) 
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Figure  4.4  Detection  probability  versus  slgnal-to-noise  voltage  ratio 
for  optimal  processing  for  various  expansion  sizes  n.  The  probability 
of  false  alarm  Is  set  to  0.05.  (taken  from  36) 


Increases  the  performance  of  the  algorithm  Improves.  As  mentioned  above 
the  performance  as  shown  Is  also  Independent  of  the  relative  strength  of 
the  clutter  noise  In  the  different  cells  but  depends  only  on  the  signal 
to  noise  ratio  In  the  target  cell.  Thus  If  we  had  only  clutter  with 
varying  parameters  along  the  A  line  the  processed  data  would  have  the 
same  statistical  properties  for  all  the  cells  containing  only  clutter. 

In  general  signal  processing  algorithms  used  to  reduce  clutter  do 
not  give  results  which  are  Independent  of  the  clutter  noise  In  the  pro¬ 
cessed  range  cell.  Thus  In  the  "Minimization"  algorithm,  analyzed  In 
Chapter  3  the  variance  of  the  processed  clutter  data  Is  a  function  of 
the  uolse  variance  of  the  cell  before  processing  and  therefore  Pf  strongly 
depends  on  the  nature  of  the  fluctuations  of  the  clutter  noise  p(h^) . 
Note  that  for  the  optimal  processing  described  by  eq.  12  Pf  Is  Independent 
of  the  fluctuations  in  the  clutter  noise  variance  p(ki)  since  In  this  case 

0»  9» 

P(f|kj^)  ■  Pj  ■  constant  so  that  Pf  -  ^  p(f |kj^)p(kj^)  ■  Pf  ^  pC^f)  ■  Pf* 

1*1  1*1 

In  figures  5(A)-(D)  the  performance  of  the  minimization  algorithm  Is 
compared  with  that  of  the  optimum  receiver  on  the  basis  of  detection 
properties  for  different  values  of  the  ratio  k  of  clutter  standard  devia¬ 
tion  Og  In  the  'worst'  cell  of  the  A  line  to  the  clutter  Of  in  the  target 
range  cell.  Each  graph  corresponds  to  a  certain  value  of  k.  We  set 
p(f|ki)  In  all  the  graphs  to  be  0.05. 

In  general  one  can  see  that  the  performance  of  the  optimum  receiver 
Is  Independent  of  k  because  the  'Student'  t  distribution  variance  Is  Inde¬ 
pendent  of  the  clutter  as  explained  above  whereas  the  performance  of  the 
Minimization  algorithm  Is  good  only  for  high  enough  Input  signal  to  noise 
ratio.  (However  the  signal  to  noise  ratio  obtained  by  the  Minimization 
algorithm  is  much  better  than  the  SNR  obtained  by  the  optimal  processing 
due  to  the  nonlinear  nature  of  the  Minimization  algorithm  as  Is  shown  In 
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Figure  4.5(A)  Detection  probability  versus  slgnal-to-noise  voltage  ratio 
for  optimal  processing,  the  original  deconvolved  signal  and  the  minimiza¬ 
tion  algorithm  for  Oq/oi"!-*  ®o/®1  ratio  between  the  stan¬ 
dard  deviation  of  a  cell  containing  only  clutter  and  the  variance  in  the 
target  range  cell. 
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Figure  4.5(B)  As  In  Flg«ire  5(A)  but  with  aQ/oi«2 
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Figure  4.5(C)  As  in  Figure  5(A)  but  with  Og/oi-A 
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chapter  3).  Figure  5  also  shows  curves  labelled  "original  signal".  For 
this  curve  the  LRT  function  used  was  that  obtained  by  simply  summing  the 
outputs  of  the  spectral  decomposition  components  r^  for  each  range  cell 
(without  division  by  the  estimated  standard  deviation).  The  use  of  the 
sum  of  r^'s  as  an  algorithm  la  attractive  due  to  the  simplicity  in  imple¬ 
mentation  and  with  the  simplifications  that  we  Introduce  later  the  optimum 
processing  for  n*®  is  equivalent  to  simply  thresholding  the  original 
signal  deconvolved  in  the  frequency  range  of  Interest  to  obtain  and  Pf . 

These  were  computed  using  the  equations 
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Figure  5(A)  shows  that,  as  was  pointed  out  earlier,  this  type  of  proces¬ 
sing  is  also  the  optimum  if  the  clutter  noise  power  is  a  constant  along 
the  A-line  (k®l).  It  is  clear  that  for  higher  ratios  of  k  the  performance 
of  such  an  algorithm  will  deteriorate,  as  is  seen  in  Figures  5(A)-(D). 

The  results  shown  in  Figure  5  (A-D)  may  be  summarized  as  follows. 
For  a  finite  number  of  windows  the  Minimization  algorithm,  (which  is  very 
powerful  in  enhancing  signal  to  noise  ratio)  is  superior  to  the  optimum 
receiver  scheme  in  detection  performance  for  large  input  signal  to  noise 
ratio.  This  is  due  to  the  fact  that  while  the  optimum  processing  is 
insensitive  to  the  variances  ratio  k  the  Minimization  performs  especially 
well  in  high  signal  to  noise  ratios  for  which  the  peak  target  echo  is 
practically  not  effected  by  the  low  noise  power  in  the  target  range  cell. 
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However  the  high  noise  variance  In  the  worst  performing  cell  Is  strongly 
reduced  by  Minimization  In  a  non  linear  fashion  that  provides  a  superior 
ROC.  As  the  ratio  k  Increases  the  apparent  advantages  of  the  optimal 
processor  over  the  Minimization  algorithm  and  the  unprocessed  signal 
become  more  evident.  Note  also  that  for  high  enough  ratio  (e.g  k^S)  the 
Minimization  algorithm  performance  Is  superior  to  that  of  the  unprocessed 
signal.  A  surprising  result  Is  the  fact  that  for  low  clutter  density 
fluctuations  (k^l  and  k»2)  Increasing  the  number  of  frequency  windows 
results  In  deterioration  In  the  ROC  of  the  Minimization  algorithm. 

EXPERIMENTAL  RESULTS 

In  this  section  we  describe  an  experimental  realization  of  the 
optimum  detector  algorithm.  Some  simplifications  are  made  that  allow 
easy  Implementation  of  the  technique  In  practical  applications. 

Recall  that  the  filters  In  the  K.L  expansion  are  matched  filters  with 
frequency  response  shape  of  slnir  T(fn~f  )/iT(f  n~f )  •  our  experimental 
realization  Gaussian  shape  filters  were  used  In  order  to  taper  the 
effect  of  the  sharp  edges  of  the  window.  This  use  of  non-matched  filters 
will  now  be  Justified  after  clarifying  some  assumptions  and  summarizing 
some  known  results. 

Up  to  this  point  we  have  not  concerned  ourselves  with  the  transmitted 
spectrum  shape.  In  ultrasound  Imaging  It  can  often  be  assumed,  due  to 
the  transducer  Impulse  response,  that  the  spectrum  shape  of  the  received 
signal  Is  roughly  Gaussian.  It  Is  shown  In  appendix  A  that  filtering 
such  Gaussian  signals  with  Gaussian  shape  filters  result  In  Gaussian  shape 
outputs.  If  the  frequency  spectrum  of  the  Input  signal  Is  much  broader 
than  that  of  the  filters  the  mean  frequency  and  the  standard  deviation  of 
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the  output  signals  will  be  practically  unchanged.  It  is  also  shown  that 
the  peak  amplitude  of  a  target  echo  will  appear  at  the  same  range  delay 
for  all  the  filters  whatever  their  center  frequency.  If  we  normalize 
each  of  the  filters  (an  operation  which  is  equivalent  to  whitening  the 
process  within  the  desired  frequency  range  and  practically  involves  as¬ 
signing  a  constant  multiplier  to  the  filters  so  that  their  output  power 
in  the  different  channels  will  be  identical),  then  the  peak  amplitudes 
appear  at  the  same  range  delay  and  have  the  same  amplitude. 

We  can  assume  that  the  clutter  noise  spectrum  in  the  input  of  each 
filter  is  much  wider  than  the  filter  bandwidth,  thus  for  the  sake  of  com¬ 
parison  between  the  matched  and  non-matched  filters  we  can  assume  the 
clutter  noise  to  be  white.  In  this  case  (see  for  example  (44))  the  de¬ 
tection  properules  of  a  matched  filter  at  a  specific  range  delay  depend 

as  is  well  known  on  the  parameter  d^rLl-”®-  which  is  the  parameter  that 

2No  2  a 

defines  the  ROC  performance  of  a  receiver.  Where  E  is  the  energy  of  the 
signal  to  which  the  filter  is  matched,  Nq  is  the  one  sided  spectral 
density  of  the  clutter  noise,  m  is  the  mean  of  the  matched  filter  output 
for  a  target  plus  noise  in  the  filter  input  and  a  is  its  standard  devia¬ 
tion.  For  the  detection  of  a  signal  with  additive  Gaussian  noise  the 
distance  is  a  function  of  the  signal  strength  m  and  the  noise  variance. 
For  mathematical  simplicity  assume  that  the  filter  is  matched  to  a  burst 
of  duration  AT  Instead  of  a  Gaussian  shaped  signal.  In  this  case 

E  -  A^AT 
2 

and 

d  -  -  /-Ai_  -  A_  ;  a2  -  N_B  ;  B  5  i  4.28 

4No  4NoB  2o  °  T 
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So  the  distance  ni/2a  for  a  filter  matched  to  a  burst  of  length  Is 
A/2/NqB  (see  for  example  59)  but  this  is  the  statistics  of  the  range  delay 
corresponding  to  the  expected  peak  target  echo  when  filtered  by  a  simple 
filter  of  bandwidth  b.  Its  mean  will  be  A  and  the  standard  deviation  will 
be  /  NqB  (assuming  of  course  that  the  peak  echo  amplitude  will  not  be 
reduced  considerably  due  to  filtering).  We  conclude  therefore  that  it 
is  possible  to  use  simple  filtering  Instead  of  matched  filtering  without 
a  substantial  loss  of  performance.  In  our  realization  the  maximum  likeli¬ 
hood  ratio  was  computed  from  the  estimated  mean  and  variance  from  the  n 
outputs  of  the  non-matched  filters  for  each  range  delay. 

In  the  experiments  the  target  echo  and  the  clutter  echo  were  created 
separately.  The  clutter  signal  consisted  of  the  echo  received  from  grainy 
metal,  whereas  the  target  signal  consists  of  the  echo  reflected  from  a 
flat  surface.  This  makes  It  easy  to  change  the  signal  to  noise  ratio  of 
the  received  echo.  The  received  signal  was  amplified  and  digitized  at 
the  RF  level  using  a  50  MhZ  digitizer  (Blomatlon).  The  information  from 
the  Blomatlon  buffer  was  read  into  a  PDF  11-23  computer.  An  FFT  was  then 
performed  and  the  spectrum  was  split  into  10  non-overlapping  Gaussian 
windows.  Instead  of  using  eq.  8  directly  we  adjusted  the  filters  gain 
to  obtain  a  white  spectrum  in  the  frequency  range  of  Interest  and  then 
applied  eq.  12. 

In  Figure  6(A)  a  typical  echo  from  a  target  plus  clutter  echo  is 
shown.  In  this  case  the  spectrum  of  the  received  signal  was  split  into 
10  channels  with  identical  bandwidths.  Though  there  was  some  correlation 
between  the  channels  the  amount  of  overlap  was  relatively  small  and  we 
assumed  the  channels  to  be  uncorrelated.  This  assumption  simplifies  the 
computations  needed  to  obtain  the  Likelihood  Ratio  Test.  In  Figure  6(B) 
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the  Likelihood  Ratio  Test  is  plotted  as  a  function  of  delay.  It  is  seen 
that  the  ratio  reaches  its  peak  at  the  target  range  delay  as  is  expected. 

For  target  enhancement  we  might  consider  using  the  LRT  instead  of  the 
original  signal.  This  will  eliminate  problems  of  incorrectly  presenting 
target  size  embedded  in  homogenous  (i.e.  uniform)  type  clutter  (such  as 
grains  for  example)  due  to  attenuation  or  shadowing.  In  these  cases  the 
local  signal  to  noise  ratio  is  the  right  measure  of  the  target  size. 


RESOLUTION  AND  RANGE  BIAS 

The  theoretical  evaluation  of  the  resolution  and  range  bias  for  both 
the  Minimization  algorithm  and  the  optimum  receiver  present  a  challange. 
Though  numerical  calculation  is  possible  it  was  found  easier  to  evaluate 
these  properties  empirically. 

We  examined  these  properties  by  simulation.  Two  target  echoes  at 
different  distances  were  summed  and  both  types  of  processing  were  applied 
to  the  signals.  These  distances  were  changed  from  half  a  wavelength  up 
to  several  range  cells  and  three  parameters  were  observed.  The  first 
was  the  ability  to  resolve  the  signals  by  eye,  the  second  was  the  effect 
of  the  processing  on  the  amplitude  of  the  processed  signal  and  the  third 
was  bias  (i.e  change)  in  the  distance  between  the  signals  due  to  proces¬ 
sing.  It  is  clear  that  signals  in  proximity  will  have  an  adverse  effect 
on  each  other.  For  the  Minimization  procedure  this  is  due  to  the  increased 
probability  of  one  echo  target  reducing  the  other  and  for  the  optimum  pro¬ 
cedure  it  is  due  to  increased  variance  caused  by  the  neighboring  target. 

It  was  found  that  the  resolution  is  phase  sensitive.  For  example 
for  about  half  a  wave  length  distance  between  the  signals,  which  means 
that  they  were  completely  unresolved  before  processing,  they  were  resolved 


for  the  optimal  procedure  after  processing  but  vere  very  noisy  for  the 
Minimization  algorithm,  as  can  be  seen  in  Figures  7(A), (B)  and  (C)  for 
both  the  suboptimum  receiver  and  the  minimization  algorithm.  However  for 
an  initial  distance  of  a  full  wave  length  the  signals  were  unresolveable 
(Figure  8).  This  behavior  was  found  to  be  cyclic,  i.e  for  an  even  number 
of  half  wavelengths  the  signals  were  not  resolveable  while  for  an  odd 
number  of  half  wave  lengths  the  signals  were  resolveable  for  optimal 
processing  and  noisy  and  meaningless  for  minimization.  For  a  distance 
larger  than  the  original  range  cell  (defined  by  the  bandwidth  of  the 
transmitted  signal)  the  signals  were  always  fully  resolveable  for  optimum 
processing  (Figures  9,10).  (However  for  minimization  the  signals  were 
resolved  only  for  a  dlstanct  of  8  wave  lengths  and  above).  It  is  impor¬ 
tant  to  note  that  these  empirical  results  were  obtained  for  a  certain 
system  and  two  processing  configurations,  namely  5  MHZ  transducer  and 
either,  4  or  10  frequency  windows. 

Further  work  is  needed  to  establish  the  general  resolution  properties 
of  the  optimum  receiver  and  the  minimization  algorithm.  However  it  seems 
that  there  is  no  significant  degradation  in  resolution  produced  by  these 
types  of  signal  processing  for  targets  positioned  by  more  than  one  range 
cell  apart. 

The  above  results  were  obtained  without  any  added  noise  and  it  is  of 
Importance  to  investigate  these  results  under  noisy  conditions.  A  clue 
to  the  behavior  under  noisy  conditions  can  be  learnt  from  the  effect  of 
the  processing  on  the  signal  amplitude  of  the  two  adjacent  signals.  One 
can  see  the  amplitude  of  the  processed  signals  as  a  function  of  separation 
It  is  seen  that  in  general  the  amplitude  of  the  processed  data  is  reduced 
substantially  if  the  separation  is  smaller  than  the  input  range  cell.  In 
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ORIGINAL  SIGNAL 


Figure  4.7(A)  Slmulatloa  of  the  echo  from  two  targets  1/2  a  wavelength 
apart . 


(B)  Processed  output  for  optimal  processing  of  the  signal  in  Figure  11(A}. 
Note  that  the  targets  are  resolved  but  that  there  Is  a  bias  In  the  dis¬ 
tance  between  them  and  that  their  amplitude  Is  substantially  reduced. 


(C)  Processed  output  for  the  minimization  algorithm  of  the  signal  In 
Figure  11(A)  Note  that  the  signals  are  not  resolved  and  that  the  ampli¬ 
tude  of  the  processed  data  Is  substantially  reduced. 
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Figure  4.8(A)  Simulation  of  Che  echo  from  two  targets  a  full  wave  length 
apart . 

(B)  Processed  output  for  optimal  processing  for  the  signal  In  Figure  12(A) . 
Note  Chat  Che  targets  are  not  resolvable  but  Chat  Che  amplitude  la  large. 

(C)  Processed  output  for  the  minimization  algorithm  for  Che  signal  In 
Figure  L2(A).  Note  Chat  Che  targets  are  not  resolvable. 
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Figure  4.9(A)  Targets  separated  by  3/2  wavelengths. 

(B)  Processed  output  for  the  Optimal  Processing.  The  targets  are  re¬ 
solvable  and  there  Is  a  bias  In  the  distance  between  them. 

(C)  Processed  output  for  the  Minimization  Algorithm.  The  signals  are 
very  noisy  and  not  resolvable. 
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Figure  4.10(A)  Targets  separated  by  4  wavelengths. 

(B)  Processed  output  for  the  Optimal  Processing.  The  signals  are  fully 
resolvable  and  the  slgnal-to-nolse  ratio  Is  high. 


(C)  Processed  output  for  the  Minimization  Algorithm.  The  output  is  noisy 
and  the  signals  are  not  resolvable. 
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Figure  4.11(A)  Targets  separated  by  8  wavelength. 

(B)  Processed  output  for  the  Optimal  Processing.  The  signals  are  fully 
resolvable  and  there  Is  no  bias  In  the  distance  between  the  signals. 

(C)  Processed  output  for  the  Minimization  Algorithm.  The  signals  are 
fully  resolvable  and  there  Is  no  bias  in  the  distance  between  the  signals 
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order  to  evaluate  the  results  qualitatively  we  have  to  compare  the  ampli¬ 
tudes  of  the  processed  data  to  the  value  of  signals  embedded  In  noise 
after  processing.  For  the  case  of  the  optimum  receiver  we  can  readily 
conclude  that  If  the  signal  to  noise  ratio  Is  higher  than  A  or  5  the 
results  would  not  change  considerably  because  the  effect  of  the  signals 
on  each  other  would  be  stronger  than  that  of  the  noise.  Minimization 
does  not  perform  as  well  for  low  signal  separation.  For  example  for 
approximately  half  a  wave  length  separation  the  signal  Is  reduced  from 
31,000  to  ~  2.4  which  for  any  reasonable  signal  to  noise  ratio  In  the 
Input  Is  negllgble.  However  for  a  full  wave  length  separation  (the  tar¬ 
gets  are  not  resolveable)  the  signal  Is  reduced  to  only  1747  while  for 
a  separation  of  a  wave  length  and  a  half  the  signal  again  Is  almost 
destroyed.  For  a  separation  of  more  than  the  original  range  cell  the 
signals  will  appear  with  a  relatively  high  signal  to  noise  ratio. 
(Figures  10,  11). 

It  Is  Important  to  note  that  while  under  some  conditions  two  adja¬ 
cent  signals  will  be  separated  by  the  processing  there  will  be  a  bias  In 
the  distance  between  the  signals  after  processing.  This  bias  was  experi¬ 
mentally  found  to  be  very  large  for  low  separation  between  the  signals. 
For  separation  larger  than  the  original  range  cell  the  bias  was  found  to 
be  negligible.  An  example  of  such  bias  can  be  seen  In  Figure  8.  One  can 
see  that  the  processed  signals  are  separated  by  much  more  than  half  a 
wavelength  which  Is  approximately  the  original  distance  between  the 
signals. 

We  can  conclude  that  for  signals  which  are  separated  by  less  than 
the  original  resolution  cell  the  results  for  both  minimization  and  the 
suboptimum  receiver  are  not  reliable  and  can  result  In  almost  a  complete 


loss  of  both  signals. 
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The  reason  for  the  surprisingly  excellent  range  resolution  perfor¬ 
mance  is  due  to  the  fact  that  in  practice  each  of  the  signals  has  several 
peaks.  The  estimated  variance  at  each  of  the  range  delays  corresponding 
to  these  peaks  due  to  the  neighboring  signal  is  clearly  phase  sensitive. 
The  variance  estimate  at  each  peak  would  be  different  and  will  depend  on 
the  exact  location  of  the  peak  in  respect  to  the  neighboring  signal.  It 
is  very  likely  therefore  that  in  at  least  one  of  the  peaks  range  delay, 
the  variance  estimate  will  be  low,  making  the  ratio  large.  This  explains 
also  the  bias  in  the  range  distance  when  the  targets  are  very  close  to 
each  other.  If  the  enhanced  peak  is  not  the  main  peak  target  echo  the 
target  range  delay  of  the  processed  data  would  seem  to  be  biased. 

SUMMARY  AND  CONCLUSIONS 

In  this  chapter  we  Introduced  a  technique  for  clutter  reduction 
which  maximizes  the  ROC  performance  in  each  range  cell.  This  technique 
has  the  property  of  a  noise  riding  threshold  as  its  performance  is  inde¬ 
pendent  of  the  local  properties  of  the  clutter  (but  is  dependent  on  the 
local  signal  to  noise  ratio).  It  was  pointed  out,  that  determining  algo¬ 
rithm  for  clutter  reduction  on  an  A  line,  requires  a  knowledge  of  the 
fluctuations  in  the  power  clutter  along  that  A-llne.  In  general  each  A 
line  would  require  a  different  processor.  It  is  therefore  more  practical 
to  use  a  detection  algorithm  which  is  optimum  for  one  cell.  This  gives  a 
likelihood  ratio  test  of  m/a  ■  Student  t,  independent  of  clutter.  We 
found  out  that  as  the  fluctuations  in  the  noise  variance  increase  so  does 
the  performance  of  the  optimal  algorithm  in  comparison  to  minimization 
and  the  deconvolved  original  signal.  (Execept  for  k-1  for  which  the 
unprocessed  signal  is  the  best).  However  in  contrast  to  Minimization  the 
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optimum  receiver  requires  a  high  dynamic  range.  Depena..ag  on  the  Input 
noise  variance  the  processed  output  can  range  from  Infinity  for  zero 
noise  power  to  almost  zero  for  signals  embedded  in  noise.  The  optimum 
receiver  algorithm  obviously  provides  the  best  detection  procedure  for 
Individual  cells.  However  as  was  pointed  out  in  chapter  3  there  are 
several  features  that  determine  the  performance  of  cn  algorithm  as  a  tool 
for  target  enhancement.  We  discussed  already  the  ROC,  dynamic  range, 
resolution  and  bias  of  the  optimum  receiver.  However  one  of  the  most 
Important  features  of  an  algorithm  is  the  Improvement  in  signal-to-noise 
ratio  for  which  the  Minimization  Is  much  superior  to  the  optimum  receiver. 
(Note  that  we  did  not  calculate  the  slgnal-to-nolse  ratio  for  the  optimum 
receiver  as  the  LRT  m/a  Is  Itself  a  measure  of  the  signal  to  noise  ratio). 
SNRE  Is  most  Important  when  additional  Information  on  the  target  Is 
available.  For  example  If  the  target  is  known  to  have  a  line  structure 
(for  an  image  composed  of  A-lines)  the  contrast  of  the  line  target  to  the 
background  noise  Is  more  important  than  ROC  of  a  single  pixie.  If  several 
pixies  along  the  A-llne  are  missing  due  to  the  processing  the  picture 
will  still  be  very  clear  to  the  eye  if  the  signal  to  background  noise  is 
high.  The  eye  will  complete  the  missing  details  and  random  Isolated 
bright  pixies  can  easily  be  removed  from  the  Image.  We  showed  that  we 
can  simplify  the  system  by  the  use  of  non-matched  filtering  Instead  of 
matched  filtering  and  exemplified  the  effectives  ass  of  the  algorithm 
by  simulation.  We  also  showed  that  surprisingly  the  range  resolution  Is 
not  significantly  deteriorated  by  the  processing.  However  a  range  bias 
Is  Introduced  for  signals  separated  by  less  than  the  original  range  cell. 
The  current  experimental  system  should  be  useful  for  correctly  presenting 
targets  embedded  In  homogenous  type  clutter  as  It  would  automatically 
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correct  for  aCteauatlon  and  shadowing*  With  some  modifications  this 
system  can  be  used  for  detection  of  targets  with  arbitrary  spectrum 
embedded  In  clutter. 


CHAPTER  5 


SUMMARY,  CONCLUSIONS  AND  FUTURE  WORK 


INRODUCTION 

As  pointed  out  in  the  introduction  the  motivation  for  this  research 
stemmed  from  the  failure  of  current  techniques  to  provide  sufficient 
clutter  and  speckle  reduction.  In  this  research  work  we  improve  these 
techniques  with  post-reception  algorithms  using  the  split  spectrum  prin¬ 
ciple  . 

At  the  outset  (chapter  2)  an  analysis  of  the  backscattered  echo  from 

random  media  is  performed.  We  first  calculated  the  average  of  the  received 

echo  from  a  regular  lattice  with  random  size  scatterers  and  found  out  that 

the  received  echo  is  proportional  to  the  scatterers  size,  scatterer 

-► 

variance,  scattering  density  and  the  function  G(r,t)  (which  is  the  Impulse 
response  of  the  system) .  Later  we  analyzed  the  backscattered  echo  from 
medium  with  random  scatterer  location  and  size,  but  uniform  average  scat¬ 
tering  density  with  results  similar  to  those  obtained  for  the  regular 
lattice  case.  Then  we  studied  the  properties  of  the  reflected  signal 
from  a  medium  with  random  scatterer  location,  random  size  and  non-uniform 
scattering  density.  We  showed  that  the  received  echo  is  strongly  depen¬ 
dent  on  the  gradient  in  scattering  density  and  calculated  the  echo  average 
for  a  step  function  in  the  scattering  density.  We  confirmed  the  theoret¬ 
ical  results  qualitatively  by  expeilment  using  sponges  to  simulate  an 


inhomogenious  medium. 
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The  results  obtained  In  chapter  2  set  the  physical  basis  to  the 
scattering  models  Introduced  In  chapter  3  and  4. 

In  chapter  3  we  analyze  the  Minimization  Algorithm,  a  technique 
that  was  first  Introduced  In  1979  by  Newhouse  and  Bllgutay  and  showed  a 
remarkable  Improvement  In  f law-to-graln  echo  suppression.  We  showed 
theoretically  that  the  slgnal-to-nolse  ratio  of  the  processed  data  is 
much  better  than  that  of  the  two  other  techniques  that  it  Is  compared 
to,  without  a  significant  loss  of  dynamic  range  (but  with  an  Inferior 
Receiving  Operating  Characteristic).  We  demonstrated  the  effectiveness 
of  the  algorithm  in  the  reduction  of  clutter.  (This  was  of  course  demon¬ 
strated  also  In  Bllgutay's  work  of  1982)  and  empirically  examined  the 
resolution  performance  of  the  algorithm.  It  was  shown  that  the  resolution 
is  not  too  seriously  adversely  affected  by  the  processing. 

In  chapter  4  we  developed  the  theoretical  basis  for  the  split  spec¬ 
trum  principle.  We  showed  that  by  this  procedure  we  transform  the  random 
process  problem  into  a  random  variable  one.  We  then  constructed  the 
optimum  receiver  and  used  a  simplified  version  of  the  this  receiver  to 
construct  an  experimental  system.  We  analyzed  the  performance  of  the 
rec?  Liieoretically  and  compared  It  to  the  performance  of  the  original 
sign  without  "special  processing".  We  also  Investigated  the  resolution 
performance  of  the  algorithm  empirically  for  signals  separated  by  less 
than  a  single  channel  range  cell  as  well  as  for  signals  separated  by  less 
than  the  original  range-cell  size.  We  show  experimentally  that  the 
resolution  Is  not  adversely  affected  by  the  processing. 
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DISCUSSION 

Though  both  of  the  techniques  Introduced  in  chapters  3  and  4  l.e.  the 
Optimal  Processing  and  the  Minimization  algorithm  were  developed  for  the 
reduction  of  clutter,  there  Is  a  significant  difference  between  them. 
While  the  optimum  receiver  scheme  allows  us  automatic  detection  of  targets 
by  Introducing  a  prespecified  threshold  that  determines  the  minimum 
slgnal-to-nolse  ratio  for  which  the  receiver  decides  that  a  target  exists 
In  the  range-cell  the  Minimization  algorithm  Is  a  technique  to  enhance 
the  contrast  of  the  target  signal  to  a  background  noise. 

The  split  spectrum  processing  Introduces  new  ways  for  signal  proces¬ 
sing.  We  demonstrate  the  versatility  of  this  algorithm  later  In  this 
chapter.  The  development  of  new  techniques  should  be  easy  due  to  the 
simplicity  In  presentation  of  the  problem  (a  set  of  finite  number  of 
random  variables,  set  to  be  Identically  distributed).  The  technique  Is 
attractive  also  due  to  the  simplicity  In  Implementation  If  simple  fil¬ 
tering  Instead  of  matched  filtering  Is  used.  In  practice  there  Is  no 
need  to  digitize  the  RF  data  and  apply  FFT  to  filter  the  signal.  The 
filtering  can  best  be  done  analogically,  which  In  most  cases  Involves 
capacitors  and  Inductors  only  (Fine  tuning  Is  needed  to  control  phase 
and  amplitude  response).  Tne  post  filtering  process  can  also  be  done  In 
hardware  (and  very  cheaply  In  most  cases). 

To  show  the  versatility  and  flexibility  of  the  use  of  the  split 
spectrum  principle  we  present  two  simple  ways  to  Improve  the  algorithms 
described  In  chapters  3  and  4. 
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PHASE  REVERSALS  AND  OVERLAPPING  FREQUENCY  WINDOWS 

If  one  examines  the  Minimization  algorithm  carefully  one  would  prob¬ 
ably  decide  that  the  principle  behind  the  algorithm  Is  the  following:  If 
there  Is  no  target  In  the  range  cell  there  Is  a  high  probability  that  the 
random  noise  phasor  Is  very  likely  to  cross  the  zero  line  If  we  observe 
Its  movement  as  we  slowly  change  the  center  frequency  of  a  frequency 
window  along  the  received  signal  band.  According  to  our  approach  we  do 
not  change  the  center  frequency  of  the  frequency  windows  continuously 
but  In  steps  and  have  n  non  overlapping  windows.  In  this  case  Instead  of 
obtaining  zero  for  noise  only  we  choose  the  phasor  which  Is  closest  to  the 
zero  line  Instead  of  the  zero  we  wished  to  "catch" .  Fortunately  we  can 
determine  If  there  would  have  been  a  zero  crossing  had  we  shifted  the 
center  frequency  of  the  windows  continuously.  Thus,  if  any  of  the  n 
random  numbers  corresponding  to  a  certain  range  delay  changes  polarity 
(or  reverses  Its  phase)  l.e.  If  there  are  both  negative  and  positive 
numbers  among  the  n  random  numbers  (one  Is  sufficient),  then  It  Is  clear 
that  there  was  a  zero  crossing.  So  we  may  decide  that  If  we  observe  a 
phase  reversal  (polarity  change)  we  will  set  the  new  processor  output  at 
zero  Instead  of  the  minimum  used  for  the  Minimization  algorithm  and 
Instead  of  the  m/a  for  the  Optimal  Processing.  It  Is  clear  that  If  the 
signal  to  noise  ratio  Is  high  enough  the  effect  of  the  new  scheme  will 
not  affect  the  target  echo  significantly.  However  for  low  slgnal-to-noise 
ratio  It  would  probably  adversely  effect  the  ROC  though  Improvement  In 
slgnal-to-noise  ratio  should  still  be  high.  (Note  that  for  very  low 
slgnal-to-nolse  ratio  the  basic  model  we  used  collapses.  The  target  echo 
decreases  and  the  assumption  that  the  target  echo  Is  much  stronger  than 
that  of  each  of  the  small  point  reflectors  that  compose  the  clutter  Is 


not  valid  any  laore).  In  Figures  1  and  2  we  show  the  same  data  shown 
earlier  in  chapter  3  and  4  but  with  the  addition  of  polarity  changes. 
The  improvement  In  signal-to-noise  ratio  is  evident  for  both  techniques. 

Another  concept  is  the  use  of  overlapping  frequency  windows.  The 
concept  was  first  Introduced  by  Newhouse  et  al  but  was  never  analyzed 
theoretically . 

Here  we  use  overlapping  windows  instead  of  the  n  orthonormal  windows 
(that  the  K.L  expansion  calls  for)  to  increase  the  probability  of  catch¬ 
ing  a  zero  crossing.  The  operation  becomes  somewhat  more  cumbersome  but 
proves  to  be  very  effective  as  is  shown  in  the  theses  of  Bllgutay  (4). 
The  addition  of  frequency  windows  should  provide  also  more  information 
to  improve  the  variance  estlmmate  in  Optimal  Processing.  In  his  thesis 
Bllgutay  investigates  empirically  the  performance  of  the  Minimization 
Algorithm  for  different  degrees  of  overlap  as  a  function  of  grain  and 
flaw  size.  It  is  pertinent  to  investigate  the  theoretical  performance 
Improvement  for  overlapping  channels  for  both  Minimization  Algorithm  and 
the  Optimal  processing.  The  derivation  Is  not  straight  forward  and 
probably  only  Chemoff  Bounds  can  be  established  for  the  problem. 

SPECKLE  REDUCTION 

In  the  Introduction  to  this  report  we  mention  speckle  reduction  as  a 
topic  of  Interest  (mainly  for  applications  In  Ultrasound  Medical  Imaging). 
Up  to  now  we  have  used  split  spectrum  processing  only  for  the  reduction 
of  clutter.  We  now  show  that  split  spectrum  processing  can  be  used  for 
speckle  suppression  (multiplicative  type  noise)  also.  Following  the  same 
approach  used  in  chapter  4  we  obtain  for  the  interval  [Ti.Tf]  the  follow¬ 
ing  hypotheses  set, 
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Figure  5.1(A)  Processed  output  for  the  signal  of  Figure  4.6(A)  for  the 
Minimization  Algorithm  without  phase  reversal  algorithm  using  10  non¬ 
overlapping  frequency  windows.  (B)  Processed  output  the  Minimization 
Algorithm  with  phase  reversal  algorithm. 
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Figure  5.2(A)  Processed  output  for  the  signal  of  Figure  4.6(A)  for  the 
Optimal  Processing  without  phase  reversal  using  10  non-overlapping  fre¬ 
quency  windows.  (B)  Processed  output  for  Optimal  Processing  with  phase 
reversal  algorithm. 
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following  the  same  procedures  and  arguments,  the  LRT  obtained  by  n  Iden¬ 
tically  distributed  Gaussian  random  variables  assuming  o]^  and  02  are  known 
Is  (see  for  example  (42)  chapter  2  pp.  107), 


5.2 


In  this  case  no  decision  can  be  made  without  further  Information  since  we 
cannot  estimate  02  1&  region  1  and  vice  versa.  However  If  we  evaluated 
the  variance  for  a  certain  region  and  want  to  determine  If  a  new  region 
belongs  to  the  same  category  l.e.  has  the  same  variance  then  a  decision 
rule  can  be  constructed.  It  Is  found  that  In  general  the  sum  of  the 
squares  constitutes  the  best  detector.  The  sum  of  the  squares  might  be  a 
good  tool  also  for  Image  formation  of  an  Ultrasonic  Image  built  by  A  scans 
If  the  grey  levels  are  calibrated  according  to  the  lowest  (or  highest) 
power  region  (this  Is  of  course  done  In  any  Ultrasound  medical  Imaging 
system  on  the  market). 

As  was  previously  mentioned  it  Is  practically  Impossible  to  build 
nonoverlapping  filters  that  will  cover  the  whole  spectrum.  Sometimes 
due  to  resolution  requirements  we  are  forced  to  use  wide  bandwidth  filters 
and  If  the  specified  number  of  windows  needed  for  speckle  reduction  Is  too 
large  to  fit  Into  the  available  transmitted  spectrum  we  use  overlapping 
windows.  In  this  case  If  the  amount  of  overlap  Is  substantial  we  are 
faced  with  a  set  of  random  variables  which  are  correlated.  For  wills  case 
the  general  solution  Is  the  following  (see  (42)) 
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rt[K-l  -  r  I  y 


where  Kq  and  are  the  covariance  matrices  of  Hq  and  respectively 
defined  as. 


E[(r  r^XHi]  - 
E[(r  rt)|Ho]  -  Kq 


For  example  assume  that  we  split  the  spectrum  Into  two  overlapping 
channels  (due  to  resolution  requirements).  Let  us  assume  that  p  Is  their 
correlation  coefficient  and  that  the  two  channels  are  assumed  to  have 
Identical  power) 


_J[  Fi(f)F2(f)df 

~rrn<£)i^-‘f 


Kl  -  E[r  r^'lHo] 


likewise 


Kq  ’  E[r  r^lHiJ 


j-1  .  J _ 

0^  2  (l-p2)  _-p  1 


The  optimum  receiver  becomes. 


120 


[1^  -  i^]  ^  (rf+rf-Zprj^r^) 

l-p2  5.9 


-  i:^](l-p2)  is  just  a  constant  multiplier.  The  receiver  should  be 
®o  ‘'1 

constructed  as  shown  in  Figure  3. 

If  we  assume  that  N  consecutive  points  on  the  A  line  are  taken  from 
the  same  statistics  but  are  independent,  the  best  detection  scheme  can 
be  shown  to  be  a  moving  average  of  the  output. 


N 

z  -  I  zi 
i-1 

which  is 


N 

z  -  j;  rJi  +  r|i  -  2p 
1«1 


5.10 


5.11 


The  first  two  terms  correspond  to  envelope  detection  and  the  next  term 
to  a  lowpass  version  of  the  correlated  signal. 


SUGGESTED  FUTURE  WORK 

In  this  section  we  outline  future  work  that  could  evolve  from  this 
research. 

1.  Use  of  coherent  reflection  from  random  media  for  the  estimation  of 
scatterer  concentration.  As  was  shown  in  chapter  2  there  are  two  ways  for 
a  coherent  reflection  to  occur,  i.e  field  gradient  and  average  scattering 
density  gradient.  The  major  problem  in  constructing  strong  enough  field 
gradient  lies  in  the  proper  design  of  a  suitable  transducer  shape  that 
will  maximize  the  field  gradient  (or  the  integral  of  as  defined  in 


chapter  2). 
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The  use  of  a  scatterlnj  concentration  gradient  to  estimate  the  scat¬ 
tering  concentration  profile  should  have  applications.  Thus  studies  on 
the  behaviour  of  the  reflection  from  smooth  boundaries  In  the  human  body 
In  relation  to  the  reflection  from  within  the  tissue  should  provide 
Information  on  the  tissue  state  If  the  scattering  model  Introduced  In 
chapter  2  can  be  assumed  (l.e.  If  the  microstructure  of  the  human  tissue 
can  be  modelled  as  a  collection  of  point  reflectors).  For  example  the 
echo  from  the  myocardium  boundary  In  relation  to  the  signal  from  within 
the  myocardium  might  Indicate  an  abnormality  of  the  myocardium  (as  the 
myocardium  fibers  are  often  modelled  as  an  ensemble  of  point  reflectors). 

As  was  pointed  out  at  the  end  of  chapter  2  this  technique  can  provide 
a  way  to  estimate  surface  grain  size  In  metal  If  the  metal  can  be  Immersed 
In  a  liquid  of  similar  characteristic  Impedance. 


ANALYSIS  OF  PHASE  REVERSALS  AND  OVERLAPPING  CHANNELS 

Phase  reversals  (or  suppression  of  polarity  changes)  and  overlapping 
channels  belong  In  some  aspects  to  the  same  category  of  signal  processing. 
Both  of  the  techniques  use  the  fact  that  a  cell  containing  clutter  only 
has  a  higher  probability  of  a  zero-crossing  than  a  cell  containing  a- 
target.  Analysis  of  the  slgnal-to-nolse  ratio  for  both  the  Minimization 
algorithm  and  Optimal  Processing  with  polarity  change  suppression  Is 
needed  In  order  to  evaluate  the  Improvement  of  these  additions  to  the 
performance  of  the  algorithms.  A  careful  examination  of  the  loss  In  ROC 
Is  also  needed  to  establish  possible  p  actlcal  Implementation  as  a  tool 
In  Image  Processing.  In  preliminary  empirical  tests  for  the  evaluation 
of  the  resolution  performance  of  this  addition  to  both  the  Minimization 
algorithm  and  the  Optimal  Processing  we  found  out  that  the  resolution  Is 


L 


123 


Improved  for  both  techniques.  In  Figures  4  to  13  we  compare  the  results 
of  "suppression  of  polarity  changes"  applied  to  both  algorithms  with  the 
results  without  "suppression".  It  is  clearly  seen  that  for  all  cases 
the  signals  are  much  better  separated  with  the  suppression  than  without 
It.  It  Is  necessary  to  establish  the  exact  Impact  that  the  "suppression" 
has  on  resolution,  signal  amplitude,  dynamic  range  and  range  bias  of 
multiple  signals 

AUTOMATIC  FLAW  DETECTION 

The  technique  In  chapter  4,  "Optimal  Processing",  can  also  be  used 
for  automatic  flaw  detection.  We  can  us'%  the  variance  estimate  to  con¬ 
struct  a  threshold  that  will  determine  a  constant  false  alarm  rate.  It 
Is  also  possible  that  flaw  characterization  In  the  presence  of  clutter 
noise  would  be  possible  If  different  type  of  flaws  would  exhibit  specific 
characteristic  spectra.  In  this  case  a  receiver  with  specific  constant 
multipliers  assigned  to  the  different  filters  would  "match"  a  specific 
flaw  type  (and  a  recognition  of  flaw  type  would  be  possible  using  some 
known  techniques  In  pattern  recognition) . 

EDGE  SHARPENING 

According  to  the  analysis  presented  In  chapter  2  sharp  gradients  In 
scattering  concentration  result  In  a  coherent  reflection  that  resembles 
a  reflection  from  a  flat  surface  or  a  large  reflector.  When  a  plane  wave 
Is  Illuminating  a  random  medium  with  a  step  In  Its  scattering  concentra¬ 
tion  we  found  out  In  chapter  2  that  the  variance  at  the  range  delay 
corresponding  to  the  Imaginary  boundary  separating  the  '■.wo  regions  of 
different  scattering  concentrations  is  approximately  half  of  the  sum  of 
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Figure  5.5  Processed  output  for  two  signals  separated  by  1/2  a  wave¬ 
length  (see  Fig.  4.7(A))  for  the  Minimization  Algorithm,  with  (B)  and 
without  (A)  phase  reversal. 
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Figure  5.7  Processed  output  for  two  signals  separated  by  a  full  wave¬ 
length  (see  Fig.  4.8(A))  for  the  Minimization  Algorithm,  with  (B)  and 
without  (A)  phase  reversal. 
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Figure  5.8  Processed  output  for  two  signals  separated  by  3/2  wavelengths 
(see  Fig.  4.9(A))  for  Optimal  Processing,  with  (B)  and  without  (A)  phase 
reversal. 
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Figure  5.9  Processed  output  for  two  signals  separated  by  3/2  wavelengths 
(see  Fig.  4.9(A))  for  the  Minimization  Algorithm,  with  (B)  and  without  (A) 
phase  reversal. 
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Figure  5.10  Processed  output  for  two  signals  separated  by  4  wavelengths 
(see  Fig.  4.10(A))  for  Optimal  Processing,  with  (B)  and  without  (A)  phase 
reversal. 
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Figure  3.11  Processed  output  for  two  signals  separated  by  4  wavelengths 
(see  Fig.  4.10(A))  for  the  Minimization  Algorithm,  with  (B)  and  without 
(A)  phase  reversal. 
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Figure  5.12  Processed  output  for  two  signals  separated  by  3  wavelengths 
(see  Fig.  4.11(A))  for  Optimal  Processing,  with  (B)  and  without  (A)  phase 
reversal . 
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Figure  5.13  Processed  output  for  two  signals  separated  by  8  wavelengths 
(see  Fig.  4.11(A))  for  the  Minimization  Algorithm,  with  (B)  and  without 
(A)  phase  reversal. 


the  variances  of  Che  two  regions.  Thus  we  can  treat  this  problem  as  a 
target  embedded  In  fluctuating  medium  where  the  standard  deviation  ratio 
between  the  cell  containing  Che  target  (the  boundary)  and  Che  cells  con¬ 
taining  clutter  only  can  take  the  value  of  either  2  or  1/2  and  the  detec¬ 
tion  performance  of  the  boundary  can  be  evaluated  using  the  type  of  ROC 
graphs  presented  In  chapter  4.  If  we  use  for  example  Optimal  Processing 
there  will  be  no  difference  between  the  processed  output  of  the  signals 
returned  from  regions  1  and  2  because  m/o  Is  Independent  of  o  as  was 
shown  In  chapter  4.  Thus  for  the  boundary  range  delay  the  processed 
output  will  depend  on  m/  a  where  m  represents  the  coherent  term  of  Che 
reflected  signal.  Use  of  Che  Minimization  Algorithm  will  result  In 
Improved  contrast  of  of  the  boundary  echo  with  respect  to  both  regions 
as  the  SNR  of  the  Minimization  Algorithm  Is  much  larger  than  that  of  any 
of  the  ocher  analyzed  techniques. 


OPTIMUM  DYNAMIC  RANGE 

One  of  Che  advantages  of  the  new  processors  Is  Chat  they  do  not 
require  a  substantial  Increase  In  dynamic  range.  Though  Minimization  is 
a  biased  estimator  of  the  signal  amplitude,  for  high  signal  to  noise 
ratios  and  not  too  large  number  of  windows  this  bias  Is  very  small  and 
can  be  considered  negligible.  Thus  there  Is  no  need  for  a  substantial 
Increase  In  amplitude  dynamic  range  of  the  processed  output .  The  same 
holds  true  for  the  Optimal  Processing  If  Che  targets  are  embedded  In  homo¬ 
genous  type  clutter  and  Che  number  of  windows  employed  Is  sufficiently 
high  to  allow  a  good  estimate  of  Che  clutter  variance. 

If  we  want  to  minimize  the  dynamic  range  requirements  for  non- 
homogenous  type  clutter  and  still  gain  high  signal  to  noise  ratio  without 
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a  need  for  a  high  number  of  frequency  windows  to  estimate  the  variance 
one  can  use  the  following  scheme:  Use  the  phase  reversal  scheme  to  decide 
If  a  target  exists  In  the  range  cell,  l.e.  If  there  was  a  phase  change  or 
not.  In  the  event  that  there  Is  no  phase  reversal  among  the  n  random 
varlbles  one  can  decide  to  present  the  original  signal  Instead  of  the 
two  algorithms  suggested  In  chapter  3  and  4.  This  scheme  will  not  require 
any  Increase  In  dynamic  range  of  the  processed  output  but  will  have  the 
clutter  reduction  properties  of  the  phase  reversal  scheme. 

If  one  wants  to  Improve  the  detection  properties  of  this  scheme  one 
way  to  do  It  Is  by  allowing  for  example  one  or  two  phase  reversals  without 
deciding  on  a  zero  (l.e.  clutter)  in  the  range  cell.  This  way  we  will  In¬ 
crease  the  probability  of  detection  with  a  similar  Increase  in  probability 
of  false  alarm. 
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APPENDIX  A 


It  can  be  shown  that  for  many  different  types  of  transmitted  signal 
and  receiver,  the  peak  of  the  echo  of  a  point  target  will  not  change  Its 
delay  with  respect  tc  the  transmit  pulse,  due  to  the  process  of  splitting 
the  received  echo  spectrum.  We  now  demonstrate  this  for  the  case  of  a 
transmitted  signal  of  Gaussian  shape 

_  t^ 

g(t)  -  e  ^  cos  oigt  (Al) 

whose  echo  from  a  point  target  Is  passed  through  n  filters  with  a  Gaussian 
transfer  function. 

It  is  clear  that  the  peak  of  g(t)  will  occur  for  t-0  (or  for  t-x 
for  a  delayed  version  of  g(t)  for  a  target  located  at  x). 

The  Fourier  Transform  of  g(t)  Is  also  Gaussian 

b2(ai-u)^)2  _t,2( 

G((d)  -  /rr b[e  2  +  e  2  ]  (A2) 

The  transfer  function  of  the  n^h  Gaussian  frequency  window  can  be  written 
as  follows 

b2((jj-(i)  )2  b2(uri-(j)  )2 

_  n  n  —  n _ n 

Hn(u))  -  /2^  bn[e  2  +  ^  2  j  (^3^ 

Multiplying  eqts.  (A2}  and  (A3),  assuming  that  the  cross  products  are 
negligible  (because  the  transfer  function  at  zero  frequency  should  be 
closely  zero) ,  we  obtain  the  output  spectrum  of  the  n^^  frequency  window 
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*2  *2 
_  '’n  (m-c4)^  -  **0  (uH-u^)^ 

G„(u)  -  G(a)).Hjj(u)  -  A  [  e  ^  +  e  ] 


where 


b*2  - 

n  a 


*  0)  b^  +  w  b^ 

uu  -  °  °  ° 

^  b2  +  bj 


and  A  Is  a  proportionality  constant  which  is  Independent  of  u). 


As  the  output  of  the  gausslan  filters  is  Gaussian  their  inverse  Fourier 
transforms  will  also  be  Gaussian  with 


g(t)  -  e  ^^n  cos  oj^t 


and  all  the  peaks  will  still  be  at  t>0< 
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APPENDIX  B 

The  Likelihood  Ratio  Test  of  the  General  Gaussian  Problem 

In  this  appendix  we  derive  the  Likelihood  Ratio  Test  for  the  dif¬ 
ferent  cases  analyzed  In  chapter  4  and  5>  We  follow  here  the  derivation 
given  In  (19). 

Let  the  set  be  jointly  Gaussian.  If  we  define 


t  - 

rl 

• 

• 

E[P]  -  A  - 

mi 

• 

• 

• 

• 

J"h_ 

and 


cov[f]  -  E((f-iii)(ifT-rtT)]  ,  K 


(Bl) 


(B2) 


Then  one  can  show  that 


Pr(i^)  -  [(2Tr)^|K|l/2j-l  exp[-i(f’f-5T)K-l (?-i) ] 


(B3) 


Under  the  first  hypothesis  we  assume  that  f  Is  specified  by  a  mean 
vector  and  a  covariance  matrix.  We  denote 

-  m^ 


(B4) 

(B5) 
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Likewise 


E[?|Ho]  -  io 

Ko  - 

The  likelihood  ratio  test  becomes 


L(r)  -  - 

PrC^l^o) 


|KQ|l/2exp[-|(i*T^T)Ki(?-i^)] 


Taking  the  logarithm  we  get 


Ln  n  +  JLlqIk, 


-  lLn| 
2 


Y 


Now  the  results  used  la  chapter  4  and  5  can  be  obtained 


(B6) 

(B7) 


(B8) 


(B9) 


1.  Eq.  4.22 

Corresponds  to  equal  covariance  matrices  with  Ki=Ko»K  and  iHq^O.  In  this 
case  the  LRT  becomes 


L(?)  -  )  - 

2  2  ^ 

■ 

which  becomes 
.  Hi 

I  Y 

Ho 


(BIO) 


(BID 


14A 


which  Is  eq.  4.22 
2.  Eq.  4.25 


Corresponds  Co  Independent  components  with  equal  variances.  R]^  becomes 
diagonal  and  Ko>kK]^,  mQ-0.  In  this  case  the  LRT  becomes 


L(?)  - 


but 


1  0 

1  0 

Kj^  ■  0^ 

Kq  -  k  a2 

• 

• 

• 

_0  1_ 

• 

_0  1_ 

(B12) 


which  Is  eq,  4.25 
3.  Eq.  5.3 

Corresponds  to  mj^-m^^^O.  The  LRT  becomes 


L(?)  -  -  l[f'^(R"l-K^l)r] 

2  2  2 


(B16) 
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and  eq.  5*3  follows 


4.  Eq.  5.2 


Readily  can  be  obtained  by  using 


Kl  -  af 

1  0 

1 

• 

Ko  ®o 

1  0 

1 

• 

1 

• 

• 

_0  1_ 

_0  1_ 

(B17) 


In  eq.  5.3. 


I 
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- 

L 


For  high  scattering  concentration  and  sharp  field  gradients  only  the 
second  term  in  the  denominator  and  numerator  are  left  and  eq.  C4  becomes 


Ln 

a  *  _ 


(C5) 


|/G(r,tj^)dv  |2 

V 

P(t2) 

• 

_  |^G(r,t2)dv  1^ 

P(ti)  _ 

At 


